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SUMMARY  AND  CONCLUSIONS 

Initial  studies  conducted  by  Battelle-Colunbus  Laboratories ,  reference 
(a), identified  potential  advantages  in  using  sore  advanced  eddy  current 
techniques  in  Naval  Air  Systens  Conmand  maintenance  inspection  applications. 

If  if o— simlnr Innn  were  made  for  a  follow-on  effort  in  which  an  existing 
Inspection  problem  would  be  addressed  using  multifrequency  approaches. 

Messrs.  Eli  Nicosia  and  Richard  Me Swain  at  the  Naval  Air  Rework  Facility 
(NAVAIREWORKFAC) ,  Pensacola,  were  contacted  by  Mr.  Richard  Deitrich  of  the 
Naval  Air  Engineering  Center,  and  a  candidate  problem  was  identified  to  be 
the  detection  and  characterization  of  corrosion  between  the  T-39  fuselage 
skin  (first  member)  and  the  airframe  (longerons  and  ribs  -  second  member). 
Present  NAVAXREWORKFAC  Pensacola  eddy  current  inspection  techniques  rely 
on  the  use  of  a  relatively  high  inspection  frequency  for  the  detection  of 
second-member  outer  surface  corrosion,  that  is  outboard  surface  corrosion. 

This  is  accomplished  by  detecting  the  presence  of  a  gap  which  is  created 
between  the  aircraft  skin/structural  member  as  a  result  of  corrosion  product 
'  expansion.  Variations  in  skin/structural  spacing  can  occur  without  the 
presence  of  corrosion,  resulting  in  the  potential  for  Incorrect  calls.  In 
addition,  second-member  inner  surface  corrosion,  that  is  inboard  surface 
corrosion,  cannot  be  detected  using  present  techniques. 

The  objective  of  the  present  effort  was  to  implement  multifrequency 
techniques  with  the  specific  goals  of: 

1.  Detecting  corrosion  on  the  inner  surface  of  the  second  member, 

2.  Quantifying  inner  surface  and  outer  surface  second-member  corrosion 
depth,  and 

3.  Differentiating  between  inner  surface  corrosion  and  the  normally 
occurring  air  gap. 

The  optimization  of  a  multifrequency  eddy  current  method  for  the 
detection  and  characterization  of  second-member  inner  and  outer  surface 
corrosion  was  modeled  using  a  computer  program  originally  developed  by 
C.  V.  Dodd  at  Oak  Ridge  National  Laboratory.  Refinements  were  made  to  the 
program  for  Increased  computational  efficiency.  The  results  of  the  analytical 
studies  showed  that  inner  surface  and  outer  surface  corrosion  could  be  detected 
and  quantified  with  regards  to  depth.  Air-gap  variations  could  also  be 
distinguished  using  multifrequency  eddy  current  data.  The  analytical  results 
were  confirmed  experimentally  using  phase-sensitive  eddy  current  design. 
Optimization  of  the  overall  eddy  current  test  design  was  greatly  aided  by 
the  use  of  the  modified  computer  program. 


Ref:  (a)  NAVAIRENGCEN  Report  NABC-92-128  of  28  Sep  1978:  Feasibility 

Evaluation  of  Advanced  Eddy  Current  Inspection  Equipment  For  Use 
in  Naval  Aviation  Maintenance  Environment 
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I.  INTRODUCTION 


A.  GENERAL. 

1.  The  objective  of  this  program  was  to  demonstrate  the  benefits  of 
multifrequency  eddy  current  technology  as  applied  to  a  particular  Naval  Air 
Systems  Command  (NAVAIRSYSCOM)  maintenance  inspection  problem — the  detection 
and  characterization  of  corrosion  hidden  beneath  the  T-39  aircraft  skin  sur¬ 
faces.  Corrosion  presently  initiates  on  the  aircraft  structural  member  at 
either  the  inner  or  outer  surface.  Inner  surface  structural  corrosion  (ISSC) 
is  presently  undetectable  using  existing  NAVAIRSYSCOM  maintenance  inspection 
procedures.  Outer  surface  structural  corrosion  (OSSC)  product,  if  severe, 
forces  a  separation  between  the  two  members.  Existing  NAVAIRSYSCOM  main¬ 
tenance  inspection  methods  rely  on  the  use  of  an  eddy  current  technique  for 
the  detection  of  this  skin/structure  air  gap.  The  gap  magnitude  is  indicative 
of  corrosion  severity.  Variation  in  skin/structure  spacing  can  occur  without 
any  structural  corrosion  being  present.  Hence  regions  which  do  not  exhibit 
any  corrosion  can  be  incorrectly  Identified  as  being  corroded.  Figure  1 
illustrates  T-39  ^SSC  and  OSSC  as  they  typically  occur. 

2.  In  order  to  systematically  Investigate  the  overall  problem  of  the 
detection  and  characterization  of  aircraft  structural  member  corrosion,  com¬ 
puter  programs  developed  by  C.  V.  Dodd  at  Oak  Ridge  National  Laboratory  were 
utilized.  These  programs  were  modified  for  increased  computational  efficiency 
and  were  used  to  study  two  important  eddy  current  test  parameters  which  are 
coil  excitation  frequency  and  coil  size.  ISSC  and  OSSC  were  modeled  as  ex¬ 
tended  areas  of  change  in  overall  skin/structural  member  thickness.  The 
effects  of  air-gap  variations  between  the  skin  and  structural  member  were 
handled  analytically  by  considering  a  three-layer  model  of  the  aircraft  skin 
and  structural  member  where  the  spacing  between  the  two  is  variable.  The 
results  of  the  computer  studies  suggest  that  a  two-frequency  multifrequency 
system  can  detect  and  characterize  ISSC  and  OSSC  as  well  as  distinguish  be¬ 
tween  air-gap  variations.  The  analytical  results  were  confirmed  experimen¬ 
tally  using  simulated  sections  of  corroded  panels.  Phase-sensitive  eddy 
current  instrumentation  incorporating  reflection  colls  was  used  for  data 
acquisition. 
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II.  REFLECTION  COIL  INSTRUMENTATION 


A.  COILS. 

1 .  Phase-sensitive  eddy  current  instrumentation  incorporating  reflec¬ 
tion  coils  was  utilized  throughout  the  course  of  this  program.  Features 
of  this  type  of  instrumentation  are  described  and  compared  with  the  more 
conventional  eddy  current  Instrumentation. 

B.  COIL  TYPES.  The  most  common  type  of  coil  used  in  conventional  eddy 
current  testing  is  the  parametric  coil.  A  parametric  coil  is  defined  as 
any  coil  in  which  the  eddy  current  generator  (drive  coil)  is  the  same  as 
the  detector  (pickup  coil).  An  alternate  type,  in  which  the  drive  and 
pickup  coils  are  physically  separated,  is  called  a  transformer  coil.  The 
voltage  detected  by  the  eddy  current  instrument  is  induced  in  the  pickup 
coil  by  transformer- like  action,  that  is  through  mutual  inductance  effects. 

1.  The  use  of  transformer  coils  offers  two  primary  advantages  over 
parametric  coils.  The  first  advantage  is  the  simplification  of  the 
associated  electronic  circuitry  afforded  by  the  isolation  of  the  power 
driving  circuits  and  the  measuring  circuits.  Secondly,  it  is  very  difficult 
to  isolate  variations  in  the  intrinsic  impedance  (R  +  ja>L)  of  parametric 
coils  from  variations  of  the  coll  impedance  due  to  the  variation  of  test 
parameters  such  as  thickness,  conductivity,  or  the  presence  of  discon¬ 
tinuities. 

2.  Variations  in  coil  resistance  with  temperature  are  especially  dif¬ 
ficult  to  eliminate,  and  may  cause  drastic  errors  in  the  eddy  current 
measurements.  Transformer  coils,  in  principle,  can  be  made  free  of  this 
effect  if  the  drive  coil  is  driven  with  a  constant  current  source,  and 
pickup  coil  voltage  is  measured  with  a  high  impedance  vector  voltmeter. 
Simplicity  of  design,  though,  is  a  considerable  advantage  of  parametric 
coils  over  transformer  coils,  as  will  be  shown  below. 

3.  The  main  contributor  to  the  design  and  use  of  transformer  coils  in 
this  country  has  been  C.  V.  Dodd  at  the  Oak  Ridge  National  Laboratory. 

His  design  of  a  transformer  coil,  which  he  calls  a  reflection  coil,  is 
shown  in  Figure  2.  As  can  be  seen,  it  consists  of  a  large  drive  coil 
surrounding  two  smaller  pickup  coils.  The  pickup  coils  are  connected  in 
series  opposition  so  that  the  voltage  induced  in  the  secondary  windings 
away  from  any  material  is  zero,  and  the  voltage  across  the  secondary  is 
due  only  to  the  electromagnetic  and  dimensional  properties  of  the  material 
being  tested. 

4.  Reflection  colls  may  be  characterized  by  several  sets  of  parameters. 
These  are: 

a.  dimension  of  the  drive  coil  (r^,  r2  »  A  1  ; ^  Parameters) 

b.  dimensions  of  the  pickup  coils  (r^»  r^,  £3.3  parameters) 

c.  location  of  the  pickup  coils  with  respect  to  the  drive  coil  (4  ^ 
1  parameter) 
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d.  electromagnetic  character  of  drive  coll  (number  of  winding  turns* 
Nj_j  1  parameter) 

e.  electromagnetic  character  of  pickup  coils  (number  of  winding  turns* 
N2;  1  parameter) 

Thus,  there  are  nine  different  parameters  that  control  the  design  of  reflec¬ 
tion  coils,  all  of  which  may  be  varied  more  or  less  independently.  Compare 
this  with  a  parametric  coil,  which  has  four  parameters  (three  dimensional 
and  one  electromagnetic),  and  it  can  be  seen  that  parametric  coils  are,  by 
far,  simpler  to  design. 

S.  Fortunately,  Dodd,  who  promulgated  the  use  of  reflection  coils,  has 
also  done  a  great  deal  to  simplify  their  design.  He  has  examined,  in  detail, 
the  interdependency  of  the  various  parameters  to  come  up  with  one  basic 
design  which  depends  most  strongly  on  only  one  parameter,  the  drive  coil 
mean  radius,  denoted  as  r  in  Figure  2.  In  general,  smaller  coils  are  used 
for  low  penetration,  high  frequency,  and  high  resolution,  and  large  coils 
are  for  deep  penetration,  low  frequency,  and  reduced  resolution. 

C.  PHASE-SENSITIVE  DETECTION.  Reflection  coils  are  typically  used  with 
phase-sensitive  detection  schemes.  Advantages  and  disadvantages  of  this 
approach  are  now  considered. 


1.  Any  sinusoid  can  be  characterised  by  three  parameters,  one  of  which 
is  frequency.  The  other  two  parameters  are  the  magnitude  of  the  signal 
and  its  phase  lead  or  lag  relative  to  some  reference  voltage.  These  will  be 
called  A  and  respectively.  These  two  parameters  may,  in  turn,  be  related 
to  two  others,  commonly  called  the  inphase  and  quadrature,  designated  x  and  y 
respectively.  The  transformation  between  these  two  sets  of  parameters  is  given 
by 


x  *  A  cos  $ 
y  ■  A  sin  $ 

tan 

A  -  (x2  +  y2)!/2 


(1) 


Detection  schemes  may  be  devised  which  measure  either  of  the  two  sets. 

These  measurements  are  usually  made  independently  of  each  other,  though, 
since  the  electronic  circuitry  necessary  for  the  above  transformation  is 
more  complex  than  the  circuitry  which  measures  either  set.  Circuits  which 
measure  are  called  phase-sensitive  instruments,  and  those  which  measure  x 
and  y  are  called  vector  voltmeters.  Measurements  of  the  magnitude  A  are  rather 
trivial,  and  will  not  be  discussed. 
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2.  Since  there  are  two  sets  of  parameters  which  describe  the  eddy  current 
response,  the  question  arises  as  to  which  should  be  used.  The  answer  to 

this  question  depends  not  only  on  the  material  parameter  being  measured,  but 
also  on  whether  the  eddy  current  system  being  used  balances  out  certain  con¬ 
stant  (or  almost  constant)  voltages  which  appear  across  the  pickup  coil. 

3.  The  voltage  across  the  pickup  can  be  separated  into  two  parts,  as 
illustrated  in  Figure  3:  the  voltage  due  to  the  mutual  inductance  between 
the  drive  and  pickup  coil  (self-inductance  in  parametric  coils)  and  the 
voltage  due  to  the  presence  of  the  test  material.  All  systems  provide  a 
means  of  nulling  or  balancing  out  the  first  of  these  two  voltages;  some 
internally  through  a  balance  network,  and  some  externally  through  the  coil 
configuration.  Typically,  parametric  coils  are  balanced  internally  and  trans¬ 
former  coils  externally,  though  this  is  not  necessarily  so. 

4.  The  remaining  voltage,  then,  is  due  to  the  presence  of  the  material 
and  is  determined  by  the  bulk  material  properties  (resistivity,  conductivity, 
and  geometry).  This,  too,  can  be  divided  into  two  parts:  one  due  to  the 
nominal  material  properties,  and  the  other  due  to  a  change  in  one  or  more  of 
these  properties.  This  is  illustrated  in  Figure  4,  where  the  sinusoid  due  to 
the  nominal  material  characteristics  is  the  vector  A,  and  that  due  to  a  change 
is  the  vector  B.  In  an  unbalanced  system,  the  two  add,  giving  a  resultant 
vector  A.  In  a  balanced  system,  A  is  eliminated,  leaving  B.  In  general, 
phase-sensitive  instruments  do  not  allow  this  kind  of  balancing,  while  vector- 
voltmeter-based  systems  do. 


5.  The  output  of  a  phase-sensitive  detector  is  a  dc  voltage  propor¬ 
tional  to  4>  ,  so  it  is  limited  to  something  corresponding  to  i  180°,  which  means 
that  0.02°  is  about  1  millivolt  on  a  ±  10-volt  full  scale  range.  Phase  measure¬ 
ments,  then,  have  a  limited  dynamic  range  regardless  of  the  amplitude  of  the 
sinusoid  being  measured.  However,  they  also  have  a  limited  sensitivity.  It 
is  not  difficult  to  show  that  the  change  in  phase  angle  ,  due  to  the 
presence  of  the  vector  B,  is  _  4 


Ap 


tan 


(2) 


where  A,  B,  6 ,  and  4$  are  defined  in  Figure  4.  A/B  ratios  of  a  thousand  or 
more  are  not  unrealistic,  which  means  that  Is  limited  to  t  0.06°,  depending 
on  0  .  On  the  other  hand,  phase-sensitive  Instruments  can  be  made  more  insen¬ 
sitive  to  liftoff  variations  than  can  vector-voltmeters,  though  a  demonstration 
of  this  is  beyond  the  scope  of  this  discussion. 


6.  The  conclusion  is  that  phase-sensitive  detectors  and  vector  voltmeters 
are  both  useful  under  different  test  situations.  Phase-sensitive  instruments 
are  good  for  measuring  bulk  properties,  surface  snd  slightly  subsurface  defects 
where  A/B  is  smaller,  and  when  sensitivity  to  liftoff  is  of  concern.  Vector - 
voltmeters,  on  the  other  hand,  are  best  for  detecting  the  smaller  volume 
flaws  and  flaws  further  into  the  material  where  A/B  is  large.  There  is,  of 
course,  no  well  defined  delineation  between  the  two,  and  there  is  much  overlap 
in  their  useful  test  applications. 
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FIGURE  3 


a)  Before  mutual  inductance  is  balanced  out 


b)  After  mutual  inductance  is  balanced  out 

VECTOR  REPRESENTATION  OF  MUTUAL  INDUCTANCE  VOLTAGE 
ACROSS  AN  EDDY  CURRENT  COIL 
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a)  In  an  unbalanced  system 


[GURE  4.  VECTOR  REPRESENTATION  OF  VOLTAGE  DUE  TO  A  TEST 
ARTICLE  NEAR  AN  EDDY  CURRENT  COIL 
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III.  EDDY  CURRENT  TEST  OPTIMIZATION  USING  COMPUTER  MODELING 
A.  INITIAL  CONDITIONS. 

1.  The  computer  program  uaed  to  model  the  airframe  corrosion  problem  is  a 
modified  version  of  C.  V.  Dodd's  program  MULTIT.  The  modifications  were  made 
to  simplify  the  data  input  routines,  and  to  decrease  the  execution  time  of 
the  integration  subroutines  by  approximately  15  percent.  The  program  is  de¬ 
signed  to  analyze  thickness  changes  in  any  given  layer  in  a  multilayer  con¬ 
ducting  medium. 

2.  As  shown  in  the  previous  section.  Dodd  has  eliminated  much  of  the 
tedium  in  designing  reflection  coils,  channeling  the  present  effort  into 
selecting  only  two  test  parameters,  the  coil  mean  radius  and  the  test  frequency. 
These  are  determined  solely  by  the  dimensional  and  electromagnetic  properties 
of  the  test  samples  under  consideration. 

3.  In  general,  there  are  two  parameters  which  characterize  any  layer  of 
the  multilayer  system  (assuming  that  all  layers  are  nonmagnetic);  these  are 
thickness  and  conductivity.  Since  the  materials  under  consideration  may  be 
divided  into  three  parts  (that  is  the  skin,  the  structural  member,  and  the 
sealant/air  gap*  or  corrosion  products  between  the  two),  there  are  six  param¬ 
eter  variations  which  should  be  investigated  to  completely  describe  the  given 
problem.  Actually  there  are  only  five,  since  the  conductivity  of  air,  sealant, 
and  corrosion  products  is  zero,  for  all  practical  purposes.  A  proper  inves¬ 
tigation  of  all  permissible  variations  was  beyond  the  scope  of  this  project; 
hence,  only  certain  variations  were  allowed,  which,  it  was  thought,  would 
indicate  the  method  of  solution  necessary  for  any  variation  of  the  above  param¬ 
eters. 

4.  The  following  assumptions  were  made: 

a.  The  thickness  of  the  skin  was  0.050  inch,  and  did  not  vary. 

b.  The  thickness  of  the  structural  member  was  0.250  inch,  and  did  not 
vary,  except  when  corrosion  was  present. 

c.  The  resistivity  of  the  skin  was  the  sane  as  that  of  the  structural 
member,  neither  of  which  varied.  This  resistivity  was  set  at  5.0 
ohm-cm,  a  nominal  value  for  aluminum. 

d.  Corrosion  of  the  structural  member  at  the  faying  surface  effec¬ 
tively  increased  the  thickness  of  the  air  or  sealant  gap  by  the 
thickness  of  the  corrosion  products,  while  decreasing  the  thick¬ 
ness  of  the  structural  meafcer  by  the  same  amount. 


*  The  terms  sealant  thickness  and  air-gap  thickness  are  used  interchangeably, 
since,  from  an  eddy  current  testing  standpoint,  they  are  equivalent. 
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e.  Corrosion  of  Che  structural  member  on  the  side  away  from  the 

faying  surface  decreased  the  thickness  of  the  structural  member, 
but  did  not  change  the  thickness  of  the  air  or  sealant  gap. 

S.  During  the  course  of  this  project ,  ten  cases  were  analyzed,  as 
detailed  in  Table  I.  OSSC,  outer  surface  corrosion,  is  taken  to  be  at  the 
structural  member  faying  surface;  and  ISSC,  inner  surface  corrosion,  refers 
to  corrosion  of  the  structural  member  away  from  the  faying  surface.  These 
ten  cases  can  be  divided  into  four  categories: 

a.  The  nominal  structure  (case  1) 

b.  Corrosion  on  the  inner  surface  of  the  structural  member  (cases  2, 
3,  and  4) 

c.  Corrosion  on  the  outer  surface  of  the  structural  member  (cases  5, 
6, and  7) 

d.  Changes  to  the  sealant  or  air  gap,  an  extraneous  variable  (cases 
8,  9,  10). 

A  pictorial  representation  of  the  four  categories  is  shown  in  Figure  5. 

TABLE  1.  GEOMETRIES  OF  THE  COMPUTER-ANALYZED  CASES 


Case 

mum 

W3MU 

39 

m 

Defect  Type 

Sealant 

ft 

iin 

1 

0.250 

0.005 

0.050 

Nominal  Structure 

2 

0.2375 

0.005 

5%  Corrosion  -  ISSC 

3 

0.225 

0.005 

0.050 

10%  Corrosion  -  ISSC 

4 

0.200 

0.005 

20%  Corrosion  -  ISSC 

5 

0.2375 

0.175 

5%  Corrosion  -  OSSC 

6 

0.225 

0.030 

I'Htvi 

10%  Corrosion  -  OSSC 

7 

0.200 

0.055 

20%  Corrosion  -  OSSC 

Sealant  Thickness 

Change  Equivalent  To: 

8 

0.250 

0.175 

IjMl  51 

5%  Corrosion 

9 

0.250 

0.030 

10%  Corrosion 

0.250 

0.055 

H 

20%  Corrosion 

6.  Existing  computer  models,  including  the  one  used  here,  have  certain  lim¬ 
itations  on  the  size  of  modeled  defects,  such  as  corrosion,  due  to  the  mathema¬ 
tical  complexity  of  modeling  electromagnetic  phenomena.  Presently,  flaw  sizes 
must  be  limited  to  either  very  large  volume  or  very  small  volume,  as  related 
to  the  size  of  the  inspection  coil.  Large  volume  analysis  was  selected  in  this 


:  j  Category  A 

Category  B 

,  Nominal  Structure 

Inner  Surface  Corrosion 

'******************* 


Category  C 

Outer  Surface  Corrosion 


Category  D 

Extraneous  Variable  (Air  Gap) 


Key 

mi  S*n 

|  1  Structural  Member 

Sealant /Air  Gap 
Corroded  Region 


FIGURE  5.  PICTORIAL  REPRESENTATION  OF  THE  VARIOUS  CORROSION 
CATEGORIES 
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project,  and  it  was  assumed  that  the  trends  predicted  by  the  large  volume 
analysis  would  accurately  follow  the  trends  of  interned late-sized  flaws.  For 
the  nost  part,  this  aasuaptlon  was  shown  to  be  true. 

B.  RESULTS  OF  THE  ANALYTIC  STUDIES.  This  section  contains  the  results  of 
the  canputer  Modeling  for  the  various  categories  detailed  in  the  previous 
section.  In  particular,  it  is  shown  how  extraneous  variables,  such  as 
changes  in  the  air-gap  thickness,  nay  be  discriminated  from  corrosion  at 
the  outer  surface  by  using  multiple  Inspection  frequencies. 

1.  DETECTION  AND  CHARACTERIZATION  OF  CORROSION  -  ISSC. 

.  a.  Cases  2,  3,  and  4  of  category  B  were  examined  to  determine  the 
relationship  between  eddy  current  phase  angle  and  test  frequency 
for  different  depths  of  corrosion.  The  results  of  the  computer 
analysis  are  shown  in  Figure  6.  Only  one  coil  was  analyzed  for 
this  case,  the  300A,  at  frequencies  ranging  from  100  Hz  to  1000 
Hz.  The  various  corrosion  depths  are  Indicated  as  a  percentage 
of  the  structural  member  thickness.  The  trends  to  be  noted  from 
these  results  are: 

(1)  As  frequency  decreases,  sensitivity  Increases,  then  decreases, 
indicating  a  point  of  maximum  sensitivity  at  about  200  Hz. 

(2)  Phase  roll-off  with  frequency  appears  to  drop*  faster  at  the 
high  frequencies  than  the  low;  that  is,  it  is  better  to  be 
below  the  maximum  than  above  it. 

(3)  The  frequency  chosen  for  this  coll  appears  to  relate  well 
with  skin  depth  6,  since  6  at  200  Hz  is  about  0.31  inch. 

b.  In  most  eddy  current  testing,  it  is  as  important  to  characterize 
flaws  as  it  is  to  detect  them;  in  this  case,  one  should  be  able 
to  give  a  figure  for  the  remaining  thickness  after  corrosion. 
Apparently  this  is  possible,  since  there  is  a  monotonically 
Increasing  relationship  between  phase  angle  and  corrosion  depth 
(see  Figure  7).  However,  it  must  be  remembered  that  this  de¬ 
rived  relationship  is  for  large  volume  flaws,  add  that  smaller 
volume  flaws  will  give  smaller  phase  angles,  at  the  same  depth, 
than  large  volume  flaws.  Thus,  any  practical  calibration  curves 
must  measure  phase  angles  for  corroded  regions  of  different 
areas  as  well  as  different  depths. 

c.  Phase  angle,  rather  than  magnitude  of  the  eddy  current  signal, 
is  used  as  the  discriminating  signal  for  two  reasons.  First,  it 
is  directly  relatable  to  the  thickness  of  the  test  piece,  and 
second,  it  can  be  made  relatively  insensitive  to  liftoff  effect, 
depending  on  the  coll  size,  frequency,  and  instrumentation  used. 

The  magnitude  of  the  eddy  current  signal,  on  the  other  hand, 

is  almost  directly  proportional  to  liftoff,  and  cannot  be 
compensated  for  liftoff  variations. 
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FIGURE  6.  DEFECT  SENSITIVITY  FACTOR  VERSUS  FREQUENCY  FOR  INNER  SURFACE  CORROSION  (300A  COIL) 
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2.  DETECTION  AMD  CHAgACTgRIIATIOH  OF  COMOSIOW  -  OggC.  Cum  5,  6  and  7 
of  category  C  (rafar  to  Figaro  5)  wore  examined  to  determine  tha  relationahlp 
between  addy  currant  phaaa  angle  and  teat  frequency  for  different  depths  of 
corroeioo  on  the  outer  aurfaca  of  tha  atructural  aatber.  Too  colls  vara  uaeft 
la  this  analysis:  tha  300A  at  frequencies  ranging  from  100  Rs  to  5  kHa,  and 
the  1S0A  at  frequencies  ranging  froe  100  Ha  to  10  kHa.  Tha  resulting  data 
for  tha  300A  la  shown  la  Figure  8*  and  that  for  the  150A  is  shown  in  Figure  9. 

Tha  treads  to  he  noted  are: 

a.  As  frequency  decreases,  the  sensitivity  Increases  then  decreases , 
pointing  out  a  region  of  uaxlaun  sensitivity  at  about  2  kHa,  for 
both  coils. 

b.  Again,  phase  rolloff  Is  higher  for  higher  frequencies,  indlca~ 
ting  that  It  la  better  to  be  too  low  in  frequency  than  too  high. 

c.  Optimal  test  frequency  correlates  well  with  skin  depth  6  since  ft 

at  2  kHa  -  0.10  inch,  just  over  the  thickness  of  the  aircraft  akin. 

Again,  it  appears  that  flaw  characterisation,  as  well  as  detection,  is  possible 
aa  long  as  flaw  area  is  taken  into  account.  Computer-generated  calibration 
curves  for  both  coils  are  shown  in  Figure  10. 

3.  DETECTION  AMD  CHARACTERIZATION  OF  SEALAHT  THICKNESS  CHARGES.  The 

pr inary  detection  uechanisa  of  OSSC  is  through  the  separation  of  the  skin  from 
noncorroded  natal  of  the  structural  member.  Obviously,  then,  any  phenonena 
which  causes  a  change  in  this  separation  will  be  interpreted  as  OSSC,  even 
though  it  nay  not  be.  Naturally  occurring  chang vs  in  the  sealant  or  air-gap 
thickness  mist  sonshow  be  discrlnlnated  agelnst.  To  do  this,  cases  8,  9,  and 
10  of  category  D  were  analysed  with  the  conputer  nodellng  progress,  using  both 
the  300A  and  150A  coils,  at  frequencies  ranging  frou  100  Hs  to  10  kHa.  The 
resulting  phase  angle  change  versus  frequency,  with  air-gap  thickness  sa  a  , 
parsneter,  are  superimposed  on  the  previous  Figures  8  and  9  in  Figures  11  and 
12.  The  laportant  points  to  note  are  listed  below: 

a.  At  high  frequency  (above  1  kHs),  there  is  essentially  no  difference 
between  equivalent  sealant  thickness  changes  and  corrosion  of  the 
outer  surface. 

b.  As  frequency  decreases,  the  phase-angle  rolloff  for  sealant  thickness 
changes  is  larger  than  for  equivalent  corrosion  of  the  outer  surface, 
thus  allowing  discrimination  between  sealent  thickness  changes  end 
corrosion. 

c.  The  rolloff  is  larger  for  the  larger  coil,  indicating  that  such 
’discrimination  could  be  more  easily  performsd  with  the  300A  than 
the  150A.  This  is  due  to  the  decreased  penetration  capabilities 
of  smaller  coils. 


21 


Corrosion  depths  are  indicated  as  a  percentage 
of  the  structural  member  thickness. 
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OUTER  SURFACE  DEEEC I  SENSITIVITY  EACTOfi  VERSUS  TEST  EfiEQUENCY  FUfi  DIFFERENT 
CORROSION  DEPTHS  <150A  COIL) 


FIGURE  11.  SENSITIVITY  FACTOR  VERSUS  FREQUENCY  FOR  AIR-GAP  THICKNESS  CHANGES 
(300A  COIL) 


SENSITIVITY  FACTOR  VERSUS  FREQUENCY  FOR  AIR-GAP  THICKNESS  CHANGES 
(150A  COIL) 
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C.  TEST  IMPLEMENTATION. 


1.  It  is  now  possible  to  envision  an  eddy  current  system  that  could  perform 
the  necessary  inspection  and  its  application.  The  system  itself  would  neces¬ 
sarily  be  a  dual-frequency  system,  where  the  test  coil  is  excited  by  two  fre¬ 
quencies  simultaneously.  The  exact  frequencies  used  would  depend  upon  the 
thickness  of  the  various  metallic  members,  and  if  the  results  presented  above 
are  generalized  to  any  structure,  the  high  frequency  is  determined  essentially 
by  the  skin  thickness,  and  the  low  frequency  is  determined  by  the  combined 

skin  and  structural  member  thickness.  For  the  case  examined  here,  with  a 
0.05- inch- thick  skin  and  0. 25-inch-thick  structural  member,  the  selected 
frequencies  are  about  2000  Hz  and  200  Hz  respectively.  Each  half  of  the  dual- 
frequency  system  would  be  identical,  except  for  frequency  selection,  and  would 
consist  of  a  phase-sensitive  Indicator  (probably  a  meter) ,  a  balance  control 
(to  set  the  Initial  phase  reading  to  zero),  and  a  liftoff  compensation  null 
(to  reduce  the  effect  of  liftoff  on  the  measured  phase).  Such  Instruments 
presently  exist  as  single-frequency  instruments;  the  modifications  necessary 
to  make  a  dual  system  are  rather  straightforward. 

2.  Three  calibration  curves  are  required.  The  first  is  an  extension  of 
Figure  7,  and  gives  the  relationship  between  phase  angle  (at  200  Hz)  and  depth 
of  corrosion  on  the  inner  surface,  as  shown  in  Figure  13.  As  indicated  above, 
the  size  of  the  corroded  area  is  important,  and  several  curves,  labeled  Al, 

A2,  A3,  etc.,  are  shown  to  emphasize  this  point.  The  second  set  of  curves. 
Figure  14,  is  an  extension  of  Figure  10,  in  the  same  manner  as  Figure  13  is  an 
extension  of  Figure  7.  Again,  various  lateral  extents  of  corrosion  should  be 
considered  in  deriving  this  calibration  curve. 

3.  The  last  set  of  curves.  Figure  15,  is  necessary  for  distinguishing 
corrosion  from  sealant  thickness  changes  and  is  derived  from  Figure  11  in  the 
following  manner.  For  the  three  corrosion  depths,  the  phase-angle  change  at 
2  kHz  is  plotted  along  the  abscissa,  and  the  phase-angle  change  at  200  Hz  is 
plotted  along  the  ordinate,  as  shown  in  Figure  15,  to  obtain  the  curve  labeled 
"outer  surface  corrosion".  Similarly,  the  phase  angle  change  for  air-gap  or 
sealant-thickness  changes  at  200  Hz  is  plotted  versus  phase-angle  change  at 

2  kHz  to  obtain  the  curve  labeled  "air-gap  change".  A  third  curve  is  drawn 
between  the  two  and  is  labeled  "decision  curve".  Testing  then  proceeds  as 
described  in  the  following  paragraph. 

4.  Call  the  low-frequency  phase  indicator  meter  A,  and  the  high-frequency 
indicator  meter  B.  The  operator  continues  testing  until  meter  A  registers  an 
indication  (either  or  both  meters  may  be  equipped  with  audible  alarms,  if 
desired).  He  then  notices  whether  meter  B  also  registers  an  indication.  If 
not,  the  Indication  is  due  to  corrosion  on  the  far  side  of  the  structural 
member,  and  the  operator  records  the  Indication.  If  meter  B  also  registers, 
then  the  operator  must  decide  whether  the  Indication  is  due  to  corrosion  of 
the  outer  surface,  or  an  air-gap  change.  This  is  done  by  locating  the  phase 
angle  at  200  Hz  versus  the  phase  angle  at  2000  Hz  of  the  Indication  on  Figure 
15.  If  the  indication  lies  above  the  decision  curve,  it  is  due  to  corrosion 
on  the  outer  surface,  and  the  operator  records  the  indication.  If  it  is  below 
the  decision  line,  it  is  due  to  an  air-gap  thickness  change,  and  the  operator 
continues  testing  until  meter  A  again  responds  to  an  indication.  This  decision 
process,  which  is  rather  simple  in  form,  is  depicted  in  Figure  16.  The  applica¬ 
bility  to  computer-based  decision  processing  is  rather  obvious. 
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*  FIGURE  13.  PHASE  SHIFT  VERSUS  CORROSION  DEPTH  FOR  VARIOUS 

'  SIZES  OF  INNER  SURFACE  CORROSION 


1 


28 


Relative  Phase  Shift  at  2kHz 
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FIGURE  14.  PHASE  SHIFT  VERSUS  CORROSION  DEPTH  FOR 
VARIOUS  SIZES  OF  OUTER  SURFACE 
CORROSION 
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Relative  Phase  Shift  at  200  Hz 


IV.  EXPERIMENTAL  VERIFICATION  OF  THE  ANALYTICAL  MODEL 


Analytical  models  are  useful  because  they  allow  the  designer  to  examine  a 
large  variety  of  test  situations  with  relative  ease  and  speed.  The  accuracy 
of  any  model  Is  suspect,  however,  because  all  models,  by  their  very  nature, 
must  include  assumptions  which  may  or  may  not  be  physically  justifiable,  and 
because  the  accuracy  of  computer  calculations  is  sometimes  low.  Dodd,  the 
developer  of  the  models  used  in  this  study,  has  verified  his  results  in  a  few 
applications,  some  simple  and  others  more  complex.  The  following  is  intended 
as  a  verification  of  the  models  in  a  somewhat  complex  test  situation. 

A.  TEST  SAMPLES.  Two  general  types  of  test  samples  were  used  in  this  study. 
The  first  type  was  to  test  the  detection  and  characterization  of  simulated  OSSC 
and  ISSC,  and  the  second  was  used  to  test  the  results  for  sealant  thickness 
changes  versus  outer  surface  corrosion.  The  first  test  sample  is  depicted  ip 
Figure  17,  and  consisted  of  a  0.25-inch  piece  of  aluminum  (7075-T6)  with 
thirteen  0.5-inch-diameter  flat-bottom  holes  milled  into  the  top  and  bottom 
surfaces  of  the  plate  to  simulate  corrosion.  As  shown,  nine  holes  were  in 
areas  away  from  edges  or  other  discontinuities.  A  row  of  fastener  holes  was 
also  present  in  the  plate  and  some  simulated  corrosion  areas  were  placed  near 
these  fastener  holes  to  examine  the  effect  of  edges  on  the  detection  of  cor¬ 
rosion.  After  manufacture,  a  skin  (0.050  inch  of  6061-TO)  was  glued  to  the 
sample.  The  second  type  of  test  sample  is  shown  in  Figure  18.  There  were 
actually  five  samples  of  this  type  made,  one  of  which  was  a  balance  sample, 
fabricated  from  a  0.25-inch- thick  piece  and  a  0.050- inch-thick  piece  of 
6061-TO.  Two  of  the  four  samples  contained  a  large  area  of  simulated  corrosion 
the  remaining  two  had  mylar  placed  between  the  two  aluminum  pieces  to  simulate 
a  sealant  or  air-gap  thickness  change  comparable  to  the  corrosion  depths  of  the 
first  two  samples. 

B.  INSTRUMENTATION. 


1.  Two  eddy-current  instruments  were  used  for  the  experimental  studies. 
The  first  was  a  Tennelec  EC  501,  which  is  based  upon  Dodd's  work.  The 
Tennelec  is  also  designed  to  use  the  reflection  coils,  with  an  operating 
frequency  which  varies  from  500  Hz  to  2  MHz,  in  1-2-5  steps.  Since  the  low- 
frequency  end  was  still  too  high  for  some  of  the  required  measurements,  the 
EC  501  was  supplemented  with  a  Super  Halec,  sold  by  Halo  Instruments.  The 
frequency  range  of  this  Instrument  is  from  10  Hz  to  9.9  kHz,  depending  on  the 
coil  size,  though  the  usable  range  appears  to  be  from  70  Hz  up.  Transformer 
colls  are  also  used  with  this  instrument,  but  are  different  from  the  reflec¬ 
tion  coils  designed  by  Dodd.  For  the  purposes  of  this  report,  it  was  assumed 
that  the  predictions  generated  foT  the  Tennelec-Dodd  reflection  colls  would 
also  hold  true,  at  least  in  principle,  for  the  Super  Halec.  That  is  to  say, 
it  was  expected  that  the  general  trends  in  the  data  would  be  similar,  though 
the  numerical  values  would  be  different. 

2.  Both  instruments  are  of  the  phase-sensitive  type;  that  is,  both 
measure  the  phase  of  the  voltage  Induced  in  the  pickup  coll,  relative  to  some 
reference  voltage.  However,  in  both  Instruments,  the  output  is  only  pro¬ 
portional  to  phase  changes,  with  the  constant  of  proportionality  being 
different  and  unknown  for  both.  Thus,  all  data,  both  in  this  section  and  the 
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previous  section,  were  normalized  to  1  by  dividing  the  output  for  a  given  fre¬ 
quency  and  test  situation  by  the  maximum  obtained  during  that  series  of  tests. 

C.  RESULTS. 

1.  A  sample  panel  of  the  type  in  Figure  17  was  used  with  both  instruments 
to  examine  the  relationship  between  phase  angle  and  frequency,  for  various 
corrosion  depths  on  the  outer  surface  of  the  structural  member.  For  both 
instruments,  the  test  coil  was  scanned  across  the  plate  until  a  maximum  phase 
change  was  noted.  The  plate  contained  corrosion  depths  deeper  than  20  percent, 
as  noted  in  Figure  17,  but  only  depths  of  5  percent,  10  percent,  and  20  percent 
were  modeled;  thus,  only  these  depths  are  shown  in  Figure  19.  Also  shown  in 
this  figure  are  the  computer-generated  plots,  for  comparison.  Note  that  the 
predicted  results  agree  very  well  with  the  experimental  data. 

2.  The  same  sample  was  used  with  the  Super  Halec  only,  since  the  Tennelec 
did  not  have  a  low  enough  frequency  response  to  examine  the  relationship 
between  phase  change  and  Inspection  frequency  for  various  depth  corrosion 
areas  on  the  inner  surface  of  the  structural  member.  The  results,  along  with 
the  computer-generated  curves,  appear  in  Figure  20.  In  general,  the  agree¬ 
ment  between  the  two  is  good,  though  the  agreement  at  the  lower  frequencies 

is  not  as  good.  This  may  be  attributed  to  two  major  factors:  firstly,  large 
volume  defects  were  analyzed,  but  saukll  volume  defects  were  made  in  the  test 
sample;  and  secondly,  the  coil  type  used  was  not  exactly  like  the  coil  type 
analyzed. 

3.  To  directly  compare  the  frequency  response  for  near-side  corrosion 
with  sealant  thickness  changes,  the  areal  extent  of  the  two  must  be  the  same. 

It  is  more  difficult  to  accurately  fabricate  small  area  sealant  thickness 
changes,  hence  the  second  type  of  sample.  Figure  18,  was  made  with  large 
area  corrosion/sealant  thickness.  The  Super  Halec  was  used  in  this  test 
since  the  lowest  frequency  of  the  Tennelec  was  again  too  high  for  this  work. 

U.  The  experimental  results  are  shown  in  Figure  21;  however,  the  pre¬ 
dicted  values  are  not  shown,  in  this  case,  because  the  numerical  agreement 
was  poor.  Note,  though,  that  the  trends  are  well  predicted. 
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Frequency,  Hz 

FIGURE  19.  EXPERIMENTAL  DATA  FOR  OUTER  SURFACE  DEFECTS 


EXPERIMENTAL  DATA  FOP  INNER  SURFACE  CORROSION 


FIGURE  21.  EXPERIMENTAL  CHARACTER  1 7AT I  ON  OF  A IP -GAP/  Al ANT  THICKNESS  CHANGES 


The  following  is  a  collection  of  the  programs  used  in  the  computer 
modeling  of  thickness  measurements  of  multilayered  structures.  The  original 
reference  for  MULTIT  by  C.  V.  Dodd,  et  al,  is  contained  In  the  report 
ORNL-4384  from  Oak  Ridge  National  Laboratory  entitled  "Some  Eddy-Current 
Problems  and  Their  Integral  Solutions." 


The  appendix  is  divided  into  three  parts:  the  first  contains  a  slightly 
modified  version  of  MULTIT,  with  the  associated  subroutines,  the  second  is  a 
program  EDTRFC  (Edit  Reflection  Coll)  which  sets  up  data  files  for  various 
reflection  coils  and  test  object  parameters  used  by  MULTIT,  and  the  third  is  a 
set  of  utility  subroutines  used  by  both  programs  for  setting  up  a  file 
structure  and  for  allowing  interaction  between  the  computer  and  the  operator. 
The  latter  two  sections,  though  written  entirely  by  staff  members  of  the 
Fabrication  and  Quality  Assurance  Section  of  Battel le  Columbus  Laboratories 
for  purposes  other  than  this  project,  are  not  considered  proprietary. 

Since  these  programs  were  designed  for  use  on  our  computer  system,  it  is 
not  advisable  for  some  other  persons  to  directly  copy  them  for  use  in  their 
(different)  computer  system,  and  this  should  be  used  for  informational 
purposes  only.  Any  person  requiring  reasonable  additional  information  should 
contact  the  staff  at  BCL  for  more  details. 
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FACE  001 


0001 

0002 

0003 

0004 

0005 

0004 

0007 

0008 

0009 

OOIO 

OOU 

0012 

0013 

0014 

0015 
001 A 
0017 
0018 


c*  * 

C..  MULT IT .FOR 
C..  WRITTEN  ST  C.  V.  POOP 
C..  ORNL--APRIL.  1973 
C.. 

C..  MODIFIED  BY  D.  T.  MAYFORD 

C..  BCL"  DECEMBER.  1978 

C. . 

C..  THIS  PROGRAM  EVALUATES  THE  SENSITIVITY  TO  A  THICKNESS 
C..  VARIATION  IN  ANY  GIVEN  LAYER  OF  A  MULT I -LAYERED 
C. .  MATERIAL 
C « • 

COMPLEX  BETA0.BETA1 .BETA. V.V97.TR. GAMMA. HUT .DRIVER 
COMPLEX  PICKUP 

REAL  L2.L3.L4.L5.L4.L7.H.N3.N4.NATL 
DIMENSION  T<  10) . U( 10) .RHO( 10) tM( 10) (BETA (10) 

DIMENSION  V<2.2) »V97<2»2) » TR<2»2> »GAHHA<3) .DRIVER! 3.5) 
DIMENSION  PXCKUPO.  S>.RL(5>.  THAG(3.5).PHASE(3»S) 

DIMENSION  SHIFT (3.5) • NUT (3*5) 

DIMENSION  C0IL<13>.  AHPCA) »  HATLU0.3).  ITITL<20>.  VERB(5> 
EQUIVALENCE  <R1 .COIL! 1 > ) * (RO.AHPC 1 > > 

•  » <  T< 1 ) . HATL ( 1 1 1 > ) • (U( 1 > rHATL ( 1 (2) ) » (RHO< 1 ) »HATL< 1 .3) > 
COMMON  X.Z.01.P1 

COMMON  /Bl/  BETA0.BETA1  1 

COMMON  /B2/  R1 .R2.R3.R4.L2.L3fL4.LS.L6.N3.N4. 

•  R6.R7.R0.R9.C4.C7.V0.G5.W.F.RS 
COMMON  /B3/  HUT. DRIVER. PICKUP. AXR1 .AIR2 
COMMON  /B4/  0AGEpX1NpX0UT.XLENpTURNSpN1A.J1p 

•  PERLAY.XLAY 

DATA  VERB/ 'COIL 7 .  'ATT  'DRFT'.  'STRT'p  'EXIT'/ 

PI  -  3 • 1415926534 
RAD  -  180. /PI 
5  JKL  -  0 

C  *  * 

C..  ENTER  COIL  PARAMETERS 
C  • « 


0019 

CALL  ATTIN(2.  1.  0.  ICUR) 

0020 

READ < 2 )  ITITL.COIL 

0021 

CALL  CLOSE ( 2  > 

0022 

RS  -  (Rl+R2>*0.5 

0023 

DO  10  1-1.9 

0024 

10 

COIL< I >  -  COIL( I >/R5 

0025 

TYPE  20.  ITITL 

0024 

20 

FORMAT (1X.20A2/) 

0027 

CALL  ATT1N(2.2.0. ICUR) 

0028 

READ < 2 )  ITITL.  N9i  MATL 

0029 

CALL  CLOSE (2) 

0030 

TYPE  20.  ITITL 

0031 

TYPE  21 

0032 

21 

FORMAT < ' SENTCR  LAYER  NO.  WITN  THICKNESS  VAR.l 

0033 

ACCEPT  22. NS 

0034 

22 

FORMAT  < 12) 

0035 

NS  -  N841 

0034 

CALL  ATTIN(2.3.0. ICUR) 
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0037  READ(2)  ITITL*  AMP 

0036  CALL  CLOSE (2) 

003V  TYPE  20.ITITL 

0040  TYPE  23 

0041  23  FORMAT < '9ENTER  OPERATING  FREG.t  ') 

0042  ACCEPT  24,  F 

0043  24  FORMAT (F10.0) 

0044  T9  -  .05 

0045  IF <N9.LT • 3.0R.N6.GE «N9 .OR • N8.E0. 1 >60  TO  1110 

0047  L7  -  L3-2.0*(L4+L5> 

0048  105  U  -  2.0*PI«F 

0049  IF(JKL.NE.O)  GO  TO  115 

0051  TYPE  110 

0052  110  FORMAT  < '  N' , 13X, ' THICK .< IN) ' ,4X, 'R<M-OHM  CN)'.3X» 

•  'M, SIGMA' , 10X, 'U' ) 

0053  115  HO  ISO  I*  1 »N9 

0054  IF  <RHO< I > .6T . 1 .0E9 )  00  TO  120 

0056  MU)  *  0.5094*UU)*F*R5*R5/RH0U> 

0057  GO  TO  125 

0058  120  M<I>  -  0.0 

0059  125  IF ( JKL.NE < 0 )  GO  TO  ISO 

0061  130  TYPE  140,1  ,  T  U  )  .RHOU  )  ,MU  >  ,U<  I ) 

0062  140  FORMAT ( '  ' , 12, 12X, 1PE12 .5,2X,E12.5>2X,E12 .5.2X.0PF6. 2) 

0063  ISO  CONTINUE 

0064  T < N9 > *  0.0 

0065  TO  -  T  <N8) 

0066  T9  «  0.05 

0067  IF < JKL.NE <0>G0  TO  205 

0069  TVAR  *  100.09T9 

c «  • 

C..  PRINT  OUT  SYSTEM  PARAMETERS 
0070  TYPE  155, N8, TVAR 

0071  155  FORMAT < '  THICKNESS  VARIATION  0F',I3»'TH  LAYER  IS 

•  ,F6.2»'X'> 

0072  TYPE  570 

0073  TYPE  160, R1 ,R2»L3 

0074  160  FORMAT < '  Rl- ' *F8 . 5, 3X, 'R2- ' ,FB.5,3X, ' DRIVER  LENGTH-' ,F8. 5) 

0075  TYPE  1 70 . R3, R4 ,L4 

0076  170  FORMAT < '  R3- ' »F8.S. 3X, 'R4- ' ,F8. 5, 3X» 'PICKUP  LENGTH-' ,F8. 5 > 

0077  TYPE  1 80 , R5 , F 

0078  180  FORMAT < '  COIL  MEAN  RADIUS-' ,F8. 5, '  INCHES' ,3X, 

•  'OPERATING  FREQUENCY-' ,1PE12. 5> 

0079  TYPE  190'LS 

0080  190  FORMAT ( '  PICKUP  RECESSED' ,F8.5> 

0081  TYPE  200*L6,L2 

0082  200  FORMAT < '  MIN  LIFT-OFF- ' »F8. 5, 3X, 'LIFT-OFF  INCREMENT-' ,F8. 5) 

C.  . 

C..  INTEGRATION  BEGINS  HERE.  THE  MUTUAL,  DRIVER,  AND 

C..  PICKUP  INDUCTANCES,  ALONG  WITH  THEIR  AIR  VALUES  ARE  CALCULATED 

C. . 

0083  205  SI  -  0.01 

0084  52  -  5.0 

0085  B1  -  0.0 

0086  B2  *  S2 
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OOB'  HO  207  !■!  «3 

0088  HO  207  J-l  ,5 

ooev  hui < i » j>  •  <o.o*  0.0 

0090  HRIVER  ( I  >  J>  -  <0.0*  0.0) 

0091  207  MCKUF-II'J)  ■  (0.0*  0.0) 

009?  A1R1  »  0.0 

0093  AIR?  •  0.0 

0094  210  11  •  (B2-BD/S1 

0095  X  -  Hl-Sl/2.0 

0094  HO  490  IBXMtll 

0097  X  •  X+Sl 

0098  TEST  •  X*L3 

0099  IF (TEST .6T.20.0)  60  TO  220 

0101  U3  ■  EXP(-TEST) 

0102  60  TO  230 

0103  220  U3  •  0.0 

0104  230  W8  ■  1.0-W3 

0105  TEST  -  X*L4 

0104  If (TEST .OT. 20.0)  00  TO  240 

0108  W4  ■  EXR(-TEST) 

0109  00  TO  250 

OHO  240  W4  -  0.0 

0111  250  W9  *  1.0-U4 

0112  TEST  ■  X*L7 

0113  IF ( TEST. OT. 20.0)  00  TO  200 

0115  W7  -  EXF-(-TEST) 

0114  60  TO  270 

0117  240  W7  -  0.0 

0110  270  WS  -  EXM-XBL5) 

0119  WO  •  1.0-W7BW4 

0120  TEST  -  X*L4 

0121  IF (TEST .OT.20.0)  00  TO  40S 

0123  U4  -  E XF'( -2.0» TEST ) 

0124  Z  •  X»R2 

0125  01  ■  R2 

C..  SUBROUTINE  BESSEL  EVALUATES  THE  INTEGRAL  OF  XBJKX) 

0124  CALL  BESSEL (VAL2) 

0127  Z  -  X*R1 

0128  01  •  R1 

0129  CALL  BESBCL(VAL1 ) 

0130  Z  -  X*R4 

0131  01  >  R4 

0132  CALL  BESSEL(VAL4> 

0133  Z  -  X*R3 

0134  01  ■  R3 

0135  CALL  BESSEL (VAL3) 

0134  S3  -  S1XVAL2-VAL1 )»( VAL4-VAL3) 

0137  84  *  S1*(VAL2-VAL1)*(VAL2-VAL1> 

0138  S5  •  S1«(VAL4-VAL3)*(VAL4-VAL3> 

0139  AS  •  WA*WS*WB*W9«W0*S3 

0140  A4  ■  W4«W8«WS*64 

0141  A7  •  W4*WS»WS*W9»W9«W0SW0*SS 

0142  IF(X.GT.SO.O)  00  TO  405 

C..  THE  LOWEST  SIONIFXCANT  LAYER «N7.  It  DCTERHIHEH 
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0144  OSUH  >0.0 

0145  N  >  N9 

0146  TEST  -  X*X*1.0E-5 

0147  280  IF<M<N).LT.TEST>  GO  TO  290 

0149  BETA<N>  -  CSORT<CMPLX<X*X*M<N>  )  >/U<N> 

0150  GO  TO  300 

0151  290  BETA<N)  >  CMPLX(X.O.O) 

0152  300  QSUM  >  0SUM9REAL < BETA<N> )*T(N)/R5 

0153  IF<0SUM.GT.20.0.0R.N.£0.1>  GO  TO  310 

0155  N  -  N-l 

0150  GO  TO  280 

0157  310  N7  •  N 

C. .  THE  MATRIX  ELEMENTS  V(N7+1 »N7)L*2  <L>1*2)  ARE 
C..  CALCULATED 

0158  V(1 *  2)  *  BETA(N+1 )-BETA<N> 

0159  V<2>2)  -  BETA(N+1 >+BETA<N) 

C..  THE  TOTAL  MATRIX  V<N9*N7>  IS  CALCULATED  BETWEEN  HERE 
C..  AND  LINE  430.  THE  GAMMA  FACTOR  IS  JUST  THE  RATIO 
C..  V<N9*N7) 1 *2/V(N9*N7 ) <2.2.  SUBROUTINE  XFORM  CALCULATES 

C..  THE  TRANSFORMATION  MATRIX  T<M+1*M)I.J  AND  SUDR0UT1HF. 

C. .  MATRIX  THEN  CALCULATES  V(K*L)IfJ 
0160  IF (N8.EQ. <N7+1 ) )  GO  TO  360 

0162  BETAO  -  BETA<N7+1 ) 

0163  N6  >  N8 

0164  IF <N8.GT . <N741 ) )  GO  TO  340 

0166  N6>  N9 

0167  340  IDX1  -  N7+2 

0168  IF<IDX1.GT.N6>  GO  TO  355 

0170  DO  350  N>1DX1*N6 

0171  BETA1  >  BETA(N> 

0172  CALL  XFORM<  N»  TR»R5tUf  T  > 

0173  JUMP  «  1 

0174  CALL  MATRIX < JUMP* V*TR) 

0175  350  BETAO  -  BETA1 

0176  355  IF(N8«LT.<N7+1))  GO  TO  410 

0176  360  097(1*2)  •  V(l,2> 

0179  V97<2* 2  >  ■  V(2,2> 

0180  IF ( N8 • EO. <N9-1 ) )  GO  TO  400 

0182  N  «  NB+2 

0183  BETAO  -  BETA<N8+1 ) 

0184  BETA1  -  BETA(N) 

0185  CALL  XF0RM<N*TR»R5*U*T> 

0186  DO  370  I>1 *2 

0187  DO  370  J>1 *  2 

0188  370  V(I*J>  -  TR(I*J) 

0189  IF <N8*EQ. (N9-2) )  DO  TO  390 

0191  BETAO  ■  BETA<N8+2> 

0192  IDX1  -  N8+3 

0193  IF< IDX1 .GT.N9)  GO  TO  390 

0195  DO  380  N«1DX1*N9 

0196  BETA1  «  DETAIN) 

0197  CALL  XF0RM(N.TR*R5tU.T> 

0198  JUMP  -  0 

0199  CALL  MATRIX(JUMP*V*TR) 
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0200 

380 

BETAO  -  BETA1 

0201 

3V0 

VV7 (lil)  »  V( 1 >2) 

0202 

VV7 (2.1)  -  V<2«2> 

0203 

400 

N  -  N8+1 

0204 

BET AO  >  DETAINS) 

0200 

BETA1  -  BETA(N) 

0206 

410 

HO  430  1*1 >3 

0207 

IF(N8.LT. (N7+1 > >00  TO  430 

020V 

AI  -  I 

0210 

T < N- 1 )  •  T0*(1.0+(AI-2.0)*TV) 

0211 

CALL  XF0RMIN.TR.R5.U.T ) 

0212 

VCX. 2)  -  W7(  1.2) 

0213 

V<2.2)  -  VV7<2.2> 

0214 

JUMP  *  1 

0210 

CALL  MATRIX! JUHP.V.TR) 

0216 

IF (N8.EQ. (N9-1 ) )  GO  TO  430 

0218 

HO  420  J-1.2 

021V 

TR< J. 1 )  «  V(Jtl) 

0220 

420 

TR( J.2>  ■  VV7< J, 1 ) 

0221 

JUMP  -  1 

0222 

CALL  MATRIX! JUMP. V.TR) 

0223 

430 

GAMMA! I)  -  V( 1 . 2>/V(2.2) 

0224 

T<N8>  ■  TO 

0220 

RL(1)  -  1.0 

0226 

TEST  -  X*L2 

0227 

IF  (TEST. GT. 20.0)  GO  TO  440 

022V 

M2  *  EXK-2.04TEST) 

0230 

GO  TO  450 

0231 

440 

U2  -  0.0 

0232 

450 

HO  470  K-2.5 

0233 

AK  -  K 

0234 

TEST  -  AK»X*L2 

0235 

IF (TEST .GT .20.0)  GO  TO  460 

0237 

RL(K)  -  H2«RL(K-1> 

0238 

GO  TO  470 

023V 

460 

RL(K>  -  0.0 

0240 

4  70 

CONTINUE 

0241 

HO  480  J-1.3 

0242 

HO  480  K*1.5 

0243 

HUT(J.K)  -  HUT  ( J.K >+GAMHA! J)*RL (K >*A5 

0244 

DRIVER ( J. K  >  -  DRIVER! J.K)fGAMNA(J)*RLCK)*A6 

0245 

480 

PICKUP <  J .  K )  -  PICKUP! J.K)4GAMMA(J)*RL(K>*A7 

0246 

485 

AIR1  •  AIR1 42 . 0*S4t ( XSL3-M6  > 

0247 

AIR2  -  AIR2+S5*!4.0*<X*L4-MV)-2.0«M7«MV*MV> 

0248 

4  VO 

CONTINUE 

024V 

B1  *  B2 

0250 

B2  ■  B2+S2 

0251 

SI  »  0.05 

0252 

IF(X.LT.V.O)  CO  TO  210 

0254 

SI  *  0.1 

0255 

IF (X.LT .29.0)  00  TO  210 

0257 

SI  *  0.2 

0258 

IF  (X.LT. 39.0)  GO  TO  210 

0260 

SI  ■  0.5 
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0261  IF < X .LT . 79 .0)  00  TO  210 

C..  THE  INTEGRATION  ENOS  HERE 

0263  IF(JKL.NE.O)  GO  TO  G80 

C.  . 

C..  NEXT.  THE  INDUCTANCES.  VOLTAGES.  AND  PHASE  SHIFTS 
C..  ARE  CALCULATED  AND  PRINTED 
0265  500  CALL  CIRCT< THAG. PHASE .GO »T1 . T2) 

0266  CALL  PHASE T(TMAG. PHASE. SHIFT .VI. SET) 

0267  TYPE  510 

0266  510  FORMAT ( '  DRIVER  RES ' »4X» ' INDUCTANCE' .AX. 

*  'NO  TURNS '.6X» 'SHUNT  CAP'. 5X» 'NOR  1H  PT'> 

0269  01  «  00*T1*T1«AIR1/U 

0270  02  *  (REAL (DRIVER (2. 3) >4AIRl >/AIRl 

0271  03  *  00*T2*T2*AIR2/U 

0272  04  *  (REAL (PICKUP! 2. 3 > ) +AIR2 ) /AIR2 

0273  TYPE  520.R6.Q1 .N3.C6.Q2 

0274  520  FORMAT!'  ' . 1PE12. 5. 2X .El2 . 5. 2X.0PFS. 1 »6X. 1PE12.S. 

*  2X.0PF9.6) 

0275  TYPE  530 

0276  530  FORMAT ( '  PICKUP  RES' *4X. ' INDUCTANCE' »4X. 'NO  TURNS'. 6X. 

*  'SHUNT  CAP'.SX.'NOR  IM  PT ' > 

0277  TYPE  520.R7.03.N4.C7.G4 

0278  TYPE  540 

0279  540  FORMAT <'  DRIVING  VOLT ' ,2X. 'SERIES  RES'.AX. 

S  'AMP  GAIN' ,6X»' INPUT  IMP') 

0280  TYPE  550.VO.R0.65.R9 

0281  550  FORMAT < '  ' .F5. 1 »9X. 1PE12.5. 2X.0PFB. 1 .6X. 1PE12.5) 

0282  TYPE  560. VI 

0283  560  FORMAT < '  DISCRIMINATOR  VOLTAGE*' . 1PE12.5) 

0284  TYPE  570 

0285  RL(1)«L6 

0286  570  FORMAT ('  ') 

0287  DO  580  1-2.5 

0288  580  RL ( I )  -  L24RL(I-1> 

0289  TYPE  590. RL 

0290  590  FORMAT!'  ' .4 (F6 . 3.8X ) .F6.3) 

0291  TYPE  570 

0292  DO  600  1-1.3 

0293  TYPE  6 1 0 . ( T MAG (I.J)fJ-1.5> 

0294  TYPE  610. (PHASE (I. J). 3- 1.5) 

0295  TYPE  610.  (SHIFT ( I . J> . J-l .5) 

0296  TYPE  570 

0297  600  CONTINUE 

0298  610  FORMAT (1X.4(1PE12<5«2X)«E12.5> 

0299  CALL  SENS(SHIFT.RAD.SEN) 

C  •  • 

C..  THE  USER  HEXT  SELECTS  ONE  OF  FIVE  POSSIDILITIES 

C. . 

0300  620  CALL  COMAND ( VERB . 5 . N5 COMMAND i  '1 

0301  GO  TO  (650.780.860.5. 1110) rN5 

C.  . 

C. .  THE  FIRST  POSSIBILITY  ALLOWS  THE  USER  TO  ALTER  THE 

C..  COIL  DESIGN  »Y  INPUTTING  THE  REQUESTED  INTEGER  DATA 

650  TYPE  660 


46  (A-8) 


I 
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0303 

0304 

0305 

0306 

0307 

0309 

0310 

0311 

0317 

0313 

0314 

0315 

0316 

0317 

0318 

0319 

0371 

0322 

0323 

0324 

0325 

0326 

0327 

0328 

0329 

0331 

0332 

0333 

0334 

0335 

0336 

0337 

0338 

0339 

0340 

0342 

0343 

0344 

0345 

0346 

0347 

0348 

0349 

0350 

0351 

0352 

0354 

0355 

0356 

0357 

0358 

0359 

0340 
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660  FORMAT < '  DRIVER  WIRE  GAGE*  TURNS*  PICKUP  WIRE  GAGE'TURNS'l 
TYPE  570 

ACCEPT  670»M1A*N2A*M3A»M4A 
670  FORMAT <41 4) 

IF ( N1 A*N2A . EQ . 0 >  GO  TO  680 

GAGE  •  N1A 

XIN  «  RlfR5 

XOUT  •  R2*R5 

XLEN  >  L3*R5 

TURNS  -  N2A 

N3  -  N2A 

N1A  «  -1 

J1  ■  1 

CALL  GAGER (R6) 

GO  TO  700 

680  IFINIA.EQ.O)  GO  TO  690 
GAGE  ■  N1A 
X1N  -  R1SR5 
XOUT  ■  R2*RS 
XLEN  -  L3XR5 
J1  -  1 

CALL  GAGER (R6) 

N3  «  TURNS 
60  TO  700 

690  IF(N2A.E0.0>  GO  TO  720 
N3  -  N2A 
TURNS  -  N2A 
XIN  -  R18R5 
XOUT  •  R2*R5 
XLEN  -  L3*R5 
J1  -  0 

CALL  GAGER (R6> 

700  TYPE  710 > TURNS* GAGE* PE RLAYfXLAY  *R6 

710  FORMAT < '  DRIVER' *F6. 1  * '  TURNS  OF  9 ' *F5. 1 > 'WIRE ' • 

•  F5. 1. '/LAYER' *F5.1* 'LAYERS' . 1PE12.5* 'OHMS') 

720  I F ( N3AFN4 A •  ED • 0 )  GO  TO  730 
GAGE  -  N3A 
TURNS  -  N4A 
N4  »  N4A 
XIN  -  R3*R5 
XOUT  -  R4*R5 
XLEN  ■  L4SRS 
J1  •  1 

CALL  GAGER (R7) 

R7  -  2.0SR7 
GO  TO  750 

730  IF (N3A*E0<0>  GO  TO  740 
GAGE  •  N3A 
XIN  »  R3*R5 
XOUT  -  R4*R5 
XLEN  -  L4SR5 
J1  -  1 

CALL  GAGER (R7) 

R7  •  2.0SR7 
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0341 

N4  «  TURNS 

0342 

80  TO  750 

0343 

740 

IF (N4A.EQ.0)  80  TO  770 

0345 

N4  -  N4A 

0344 

TURNS  -  N4A 

0347 

XI N  -  R39R5 

0348 

XOUT  «  R4*R5 

0349 

XLEN  -  L49R5 

0370 

J1  -  0 

0371 

CALC  8A6ER(R7) 

0372 

R7  »  2.04R7 

0373 

750 

TYRE  740*  TURNS* GAGE  * PERLAY * XLAY  *R7 

0374 

740 

< 

FORMAT  < '  F'ICKUF"  *F4. 1  *  '  TURNS  EA9' *F5. 1 • 'WIRE' > 

I  F5.1* '/LAYER '.F5.1* 'LAYERS '.1PE12. 5* 'OHMS' ) 

0375 

770 

r* 

80  TO  500 

W  4  4 

C  4  4 

THE  SECOND  POSSIBILITY  ALLOWS  THE  USER  TO  ALTER  THE 

C  4  « 

ATTENUATOR  DESIGN  BY  INPUTTING  THE  DATA  REQUESTED  IN 

C.. 

E  FIELD  FORMAT 

0374 

780 

TYPE  790 

0377 

790 

1 

FORMAT  < '  DRIVER  SERIES  RES'  SHUNT  CAP'  PICKUP  SHUNT  RES' 
t  SHUNT  CAP ' ) 

0378 

TYPE  570 

0379 

ACCEPT  800fPlfP2'P3'P4 

0380 

800 

FORMA  T(4F10«0) 

0381 

IF (PI . EG . 0. 0)  80  TO  610 

0383 

RO  *  PI 

0384 

810 

IF(P2.EG.0.0>  80  TO  820 

0384 

C4  -  P2 

0387 

820 

1F(P3.EO.O.O)  60  TO  830 

0389 

R9  -  P3 

0390 

830 

IF(P4.EQ.0.0>  80  TO  840 

0392 

C7  «  P4 

0393 

840 

Q5  -  1 . 0/ (W*SQRT (01 *Q2*C4  > ) 

0394 

04  -  SORT < 01*02/04) 

0395 

07  -  R0/(W*Q1*Q2) 

0394 

08  -  1 • 0/ ( WtSQRT (03*04*C7) ) 

0397 

09  -  SORT (Q3S04/C7) 

0398 

010  -  R9/(W*Q3*04) 

0399 

TYPE  850 »05 *04 *07 

0400 

850 

i 

FORMAT ( '  DVR  CKT ' * 1PE12. 5* '  BELOW  RES' . 1PE12.5* 

1  '  OPT  RES''1PE12.5*'  RES/REACT.') 

0401 

TYPE  8SS'06'09'Q10 

0402 

855 

FORMAT < '  P-U  CKT' . 1PE12.5* '  BELOW  RES'  * 1PE12.5* 

»  '  OPT  RES' .1PE12.5*'  RES/REACT.') 

0403 

r* 

GO  TO  500 

L*  9  4 

c. . 

THE  THIRD  POSSIBILITY  ALLOWS  THE  USER  TO  EXAMINE  THE 

c. . 

EFFECT  OF  DRIFTS  ON  THE  CIRCUIT 

0404 

840 

TYPE  870 

0405 

870 

FORMAT ('  SYSTEM  DRIFT  VARIATIONS') 

0404 

DR1  -  0.01 

0407 

DR2  •  0.01 

0408 

DR3  -  0.01 
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040V 

DR4  »  0.01 

0410 

DRS  -  0.01 

0411 

DR4  -  0.01 

0412 

DR7  »  0.01 

0413 

DR8  -  0.01 

0414 

DR9  -  0.01 

0415 

GO  TO  890 

0414 

860 

CALL  CIRCTKSET.Vlr6HIFTtRAb.SEN> 

0417 

GO  TO  < 1 070 • 1 1 00  > • JKL 

0416 

890 

TYPE  900 

0419 

900 

FORMAT!'  XVARI AT I ON '.3*.' PARAMETER  VAR'.i: 

«  'RADIANS' »7X» 'DEGREES' »7X. 'Z  OF  RANGE 

0420 

IF(DRl.EO.O.O)  GO  TO  930 

0422 

R4  -  R4*(  1 .0+DR1 ) 

0423 

DR100  -  100.0SDR1 

0424 

TYPE  920. DR 100 

0425 

920 

F0RMAT('«'.F4.1.10X. 'DRIVER  RES  ') 

0424 

CALL  CIRCTK SET. VI .SHIFT .RAD. SEN) 

0427 

R4  -  R4/ ( 1 . 0+DR1 > 

0426 

930 

IF (DR2.E0 • 0 . >  GO  TO  950 

0430 

R7  ■  R7»(1.04DR2) 

0431 

DR100  -  100.SDR2 

0432 

TYPE  940.DR100 

0433 

940 

FORMAT < ' 4 ' >F4. 1 • 10X. 'PICKUP  RES  '> 

0434 

CALL  C1RCT1 (SET. VI .SHIFT. RAD. SEN) 

0435 

R7  -  R7/U.0+DR2) 

0434 

950 

IF(DR3.EQ.0.0)  GO  TO  970 

0436 

C4  -  C4SO.O+DR3) 

0439 

DR 100  -  100.06DR3 

0440 

TYPE  940.DR100 

0441 

940 

FORMAT ('t'.F4<1.10X. 'DVR  SHUNT  CAP') 

0442 

CALL  CIRCT1 (SET. VI. SHIFT. RAD. SEN) 

0443 

C4  -  C4/( 1 .04DR3) 

0444 

970 

IF ( DR4 . EO . 0 . 0 )  00  TO  990 

0444 

C7  -  C7*( I.04DR4) 

0447 

DR100  -  100.SDR4 

0446 

TYPE  980.DR100 

0449 

980 

FORMAT ( '4'. F4. 1.1 OX. 'P-U  SHUNT  CAP') 

0450 

CALL  CIRCT1 (SET. VI .SHIFT. RAD. SEN) 

0451 

C7  »  C7/( 1 .0+PR4 ) 

0452 

990 

IF(DR5.EQ.O . 0)  GO  TO  1010 

0454 

RO  -  ROS(l.OFDRS) 

0455 

DR100  >  100. SDRS 

0454 

TYPE  1 000 . DR 1 00 

0457 

1000 

FORMAT ( '• ' .F4 • 1 . 10X. 'SERIES  RES  ') 

0456 

CALL  CIRCT1 (SET. VI .SHIFT. RAD. SEN) 

0459 

RO  -  R0/< 1 .0+DR3) 

0440 

1010 

IF (DR4.E0.0.0)  00  TO  1030 

044? 

R9  -  R9* ( 1 . 0+DR4 ) 

0443 

DR100  »  100.6DR4 

0444 

TYPE  1020 .DR 100 

0445 

1020 

FORMAT ('t'.F4.1«10X. 'AMP  INPUT  RES') 

0444 

CALL  CIRCT1 (SET. VI .SHIFT .RAD. SEN) 

0447 

R9  »  R9/( 1 .04DR4) 

PAGE  009 
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RAGE  010 


0468 

0470 

0471 

0472 

0473 

0474 

0475 

0476 

0478 

0479 

0480 

0481 

0482 

0483 

0484 

0485 

0486 

0488 

0489 

0490 

0491 

0492 

0493 

0494 

0495 

0496 


0497 

0498 


1030  IF ( DR  7 . EQ . 0 . 0)  60  TO  1050 
VO  -  V0*< 1 . 0+DR7) 

DR 100  -  100.*DR7 
TYPE  1040.DR100 

1040  F0RMAT('*'.F4.1.10X. 'APPLIED  VOLT  ') 

CALL  CIRCT1 <  SET  r VI 'SHIFT  . RAD • SEN ) 

VO  •  V0/< 1 .0+DR7 ) 

1050  I F  <  DR8 . EQ • 0  <  0  >  GO  TO  1080 
F  -  FX1.0+DR8) 

DR 100  -  100.8DR8 
TYPE  1060 i DR 100 

1060  FORMAT <'*'»F4«ltlOX» 'FREQUENCY  '> 

JKL  -  1 
00  TO  105 

1070  F  -  F/( 1 • 0+DR8 > 

JKL  ■  0 

1080  IF<DR9.E0.0>  GO  TO  620 
RS  -  R5*< 1 . 0+DR9) 

DR100  -  100.CDR9 
TYPE  1090. DR 100 

1090  FORMAT ('4'>F4.1.10Xr 'MEAN  RADIUS  ') 

JKL  -  2 
60  TO  105 

1100  RS  -  R5/U.0+DR9) 

JKL  «  0 
60  TO  620 

C.  . 

C..  THE  FOURTH  P0SSI8ILITY  RETURNS  TO  PROGRAM  START 

C.  . 

C..  THE  FIFTH  P0SS18JLTY  ENDS  ALL  CALCULATIONS 

C.  . 

1110  CALL  EXIT 
END 
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0001  SUBROUTINE  BESSEL  <  UAL ) 

0002  COMMON  X.  Z*  01 t  PI 

0O03  IF (2 . GT .5.0)  GO  TO  1510 

0005  L  '  2.0*Z+3.0 

OOOo  FI  -  0.5*01*01*01 

0007  UAL  ■  F 1/3.0 

0008  BO  1500  I  ■  It  L 

0009  A I  *  I 

OOIO  FI  «=  -FI  *0. 25*2*2/ <AI*A1 +AI  > 

0011  1500  UAL  «  VAL+Fl/<  2.0*Al  +  3.0) 

0012  GO  TO  1520 

0013  1510  Z1  «  l./Z 

0014  XO  «  <<<-188.1357*Zl+109.1142)*Zl-23. 793331*21+2.050931 >*Z1 

0015  XO  *  ( < XO-0 . 1 730503 ) *Z1 +0 • 7034845 >*Z 1-0. 0641 09E -3 

001*  XI  •  <(<-5. 817517*Zl+2. 1058741*21-. 6896196>*Z1+ . 4952024 )*Zl 

0017  XI  =  <Xl-0.l87344E-2)*Zl+0. 7979095 

0018  UAL  ■  <1 .0-SGRT<Z)*<Xl*COS'it -PI/4. >-X0*SlN(Z-RI/4. ) > >/<X*X*X) 

0019  1520  RETURN 

0020  END 
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OOOl  SUBROUTINE  XFORM<N,  TRf  RS*  U>  T> 

0002  COrtHON/Bl /BET AO »  BETA1 

0003  COMPLEX  BETAO,  BETA1,  EX,  TR 

0004  DIMENSION  EX<2>,  U<10>.  T<10>,  TR(2,  2) 

0005  EX( 1 >  «  CEXP(-BETA0*U(N-1 >*T<N-1 >/RS> 

0006  EX (2)  -  l.O/EXU) 

0007  00  1610  I  -  1,  2 

0008  00  1610  J-l,  2 

0009  K  -  1+J 

0010  JF(K.£0,3)  00  TO  1600 

0012  TR<  1 ,  J)  -  ( BETA1 4-BET  AO)  *EX  <  J) 

0013  GO  TO  1610 

0014  1600  TR( 1 «  J)  -  (BETAl -BETAO ) *EX ( J) 

0015  1610  CONTINUE 

0016  RETURN 

0017  END 
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0001 

SUBROUTINE  MATRIX (JUMP.  V.  TR) 

0002 

COMPLEX 

V.  0.  TR 

0003 

DIMENSION  0(2.  2).  V<2.  2>t  TR<2.  2> 

0004 

IF (JUMP. 

EQ.l)  GO  TO  1720 

0006 

DO  1700 

1*1.  2 

0007 

1700 

0(1.  1) 

-  V(I.  1> 

0008 

DO  1710 

I»1 .  2 

0009 

V(I.  1) 

■  (0.0.  0.0) 

0010 

DO  1710 

J-l.  2 

0011 

1710 

V(I.  1) 

-  V(I.  1 )+TR( I ,  J)*Q(J.  1) 

0012 

1720 

DO  1730 

1-1.  2 

0013 

1730 

Q( 1 .  2) 

-  V(I.  2) 

0014 

DO  1740 

1-1.  2 

0015 

V(I»  2) 

-  (0.0.  0.0) 

0016 

DO  1740 

J-l.  2 

0017 

1740 

V(I,  2) 

-  V(I.  2)+TR( I .  J)*G(Jr  2) 

0018 

RETURN 

0019 

END 
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0001  SUBROUTINE  CJRCT (THAG.  PHASE .  00.  Tl«  T2) 

0002  COHPLEX  HUT .  DRIVER.  PICKUP.  21.  22.  23.  24.  25.  26.  27 

0003  COHPLEX  DENOH.  TNUH.  VOLT 

0004  COHHON  /B2/  Rl.  R2.  R3.  R4.  L2.  L3.  L4.  LS.  L6.  N3»  N4. 

•  K6«  R7 .  RO.  R9 .  C6.  C7.  VO.  G5.  U»  F»  RS 

0005  COHHON  /B3/  HUT.  DRIVER.  PICKUP.  AIR1.  AIR2 

0006  D I HENS I ON  HUT(3.  5).  DRIVER (3.  5).  PICKUP<3.  5) 

0007  DIHENSION  THAG< 3.  5>.  PHASE<3.  S> 

0008  REAL  L3.  L4.  N3»  N4>  L2.  L5.  L6 

0009  T 1  ■  N3/< (R2-R1 )*L3> 

0010  T2  »  N4/<(R4-R3)*L4> 

0011  QO  -  6.300475204E-7*F*R5 

0012  21  -  CHPLX ( W*C6*R0 .  -1.0) 

0013  22  -  CHPLX <W*C7*R9.  -1,0) 

0014  23  -  CHPLX (0,0.  -RO) 

0015  24  -  CHPLX <0,0.  -R9) 

0016  DO  1800  I-l.  3 

0017  DO  1800  J-l.  5 

0018  25  -  Q0*T1*T2*HUT(I»  J) 

0019  26  -  Q0*T1 *T 1 * (0 . 0 .  1 .0 >t (DRIVER ( I .  J14AIR1 > 

0020  27  -  Q0*T2*T2*<0.0»  1 , 0)* (PICKUP < 1 .  JHAIR2) 

0021  DENOH  -  21*22*25*25+<21*<26+R6)423>*<22*<27+R7>424> 

0022  TNUH  -  V0*R9*G5*<0.0»  -1.0 >*25 

0023  VOLT  -  TNUH/ DENOH 

0024  THAG(1.  J)-  CABS (VOLT) 

0025  1800  PHASEU.  J)  -  ATAN2(AIHAG(  VOLT  )  »  REAL  (VOLT  >  > 

0026  RETURN 

0027  END 
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0001  SUBROUTINE  PHASE T <  Ai  B»  C*  VI  *  03) 

0002  DIMENSION  A<3.  5).  B(3»  5).  C<3.  5) 

0003  01  -  A<2.  5  >  PSIN( B ( 2 '  5>)-A(2t  1)*SIN<B<2>  11) 

0004  02  -  -<A<2«  5>*C0S(P(2»  5)>-A(2.  1)*C0S(B<2»  1)1) 

0005  03  -  ATAN2(01 t  02) 

0006  VI  «  A  ( 2 1  1 >*SIN<03+B<2»  1)1 

0007  HO  1900  I-l*  3 

0008  HO  1900  J-l.  5 

0009  1900  C(I*  J)  ■  03-ATAN21V1 •  SQRT(A<I>  J)*A(I»  J)-V16Vl))+B(If  J> 

0010  RETURN 

0011  END 
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0001  SUBROUTINE  8£NS(TH.  RAO.  TNI > 

0002  DIMENSION  TH<3.  S) 

0003  TH9  •  THU.  1) 

0004  THO  ■  TH< 1 .  1> 

OOOS  DO  2010  1*2.  5 

0006  IF<  TH< 1 .  I )  .LT. TH9>  60  TO  2000 

0008  TM9  -  TH< 1  *  I> 

0009  2000  IF<TH< 1 »  D.GT.TMO)  GO  TO  2010 

0011  THO  -  TH< 1 «  I> 

0012  2010  CONTINUE 

0013  TH1  »  <TH< 1 .  1 >+TH( 1 .  2>+TH(l.  3>+TH(l.  4>*TH<1.  5)1/5. 0 

0014  TH1  -  TH1-(TM(3»  1>+TH<3.  2>*TH<3.  3)+TH<3.  4)+TH<3.  5))/5. 

0015  XLO  -  TH9-TH0 

0014  PERCT  -  100. tXLO/THl 

0017  DEGR1  -  TH18RAD 

0018  DEGR2  -  XLOtRAD 

0019  TYPE  2020.  TNI .  XLO.  PERCT 

0020  2020  FORMAT < '  PHASE  SHIFT*' .  1PE12.5.  2X.  'LIFTOFF*' 

•  •  E12.5.  2X»  'X-'.  E12.5) 

0021  TYRE  2030.  DEGR1.  DE6R2 

0022  2030  FORMAT < '  DEGREE! ' .  1PE13.S.  5X.  E13.5) 

0023  RETURN 


0024  END 
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FORTRAN  IV  V01C-03G  FRI  19-0CT-79  09558155  RAGE  001 

0001  SUBROUTINE  GAGER (RES) 

0002  COMMON/ 84 /G t  PI.  02.  RLN.  T5.  NIAt  Jit  0.  E 

0003  F'l  »  3.1415926536 

0004  IF(Jl.EQ.l)  GO  TO  2120 

0006  N1A  -  -1 

0007  03  >  0.954S0RT < (D2-P1 XRLN/T5) 

0008  X2  -  l«0371E-5/( 8341*3 > 

0009  0  -  40.0 

0010  2100  G  «  40. +10. *<AL0G<X2)-AL0G(. 99B9+. 017*<G/10.-1 . > ) >/2. 30259 

0011  IF<ABS(G)-G> . LT. 1 .E-4)  GO  TO  2110 

0013  0  •  G 

0014  GO  TO  2100 

0015  2110  16  -  G 

0016  G  •  IG 

0017  2120X2-  < .9989+.017*(G/10.-l. > >*10. **<G/10.-4 . > 

0018  03  -  SORT < 1 . 0371E-S/X2) 

0019  1F(G«GT .40. )  GO  TO  2130 

0021  X3  -  < . 4 606554 ALOG ( 03*1 . E3>- . 43444)61 . E-3 

0022  GO  TO  2140 

0023  2130  X3  -  ( 98 . 02228*03+2 . 56791E-2 ) *1 . E-3 

0024  2140  IP  >  <RLN/(03+X3> ) 

0025  0  -  ID 

0026  IE  -  <02-01 )/<D3+X3> 

0027  E  -  IE 

0028  IF (N1A.EQ. -1 )  GO  TO  2150 

0030  T5  -  D*E 

0031  2150  RES  -  T5*X2* ( D2+D1 )*RI /12 . 0 

0032  RETURN 

0033  ENO 
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FORTRAN  IV 


V01C-03G  FRI  19-OCT-79  09S59S22 


RAGE  001 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0010 

0011 

0012 

0013 

0014 

0015 

0016 

0017 


SUBROUTINE  CIRCT1 < TH1 •  VI.  TH.  RAG.  SEN) 

PI MENS I ON  TH< 3»  5>.  TMAG(3.  5).  PHASE <3.  5) 
CALL  C1RCT (TMAG*  PHASE.  00.  Tl.  T2> 

01  ■  0.0 
BO  2200  1*1.  3 
BO  2200  J«l.  5 

02  -  TH1-ATAN2< VI .  SORT (TNAG(I.  J>**2-V1*S2) ) 
•  +PHASE ( I »  J) -TH< I »  J> 

IP (ASS (01 ) .GT .APS (02) >  GO  TO  2200 
01  -  02 
2200  CONTINUE 

G2NEW  ■  RAB*01 
Q3  «  100.SG1/SEN 
TYPE  2210.  01.  02NEW.  03 
2210  FORMAT! IX.  2<1PE12.5.  2X>.  1PE12.5) 
RETURN 
ENB 
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FORTRAN  IV  V01C-03C  FRI  19-OCT-79  09J59S44 


FACE  001 


0001 

0002 

0003 

OOOS 

0006 

0007 

0008 

0009 

0010 

0011 

0012 

0013 

0014 

001S 

0016 

0017 

0018 

0019 

0020 


SUBROUTINE  BE6S1<VAL> 

COMMON  X.  2,  Ql>  FI 

IF  12 .GT .5.0)  CO  TO  1510 

L  -  2.0*243.0 

FI  -  0.5*2*01*01 

VAL  »  Fl/3.0 

CO  1500  I  ■  1.  L 

AI  -  I 

FI  -  -F 1*0. 25*2*2/ <AI*AI4AI ) 

1500  VAL  ■  VAL4F1/(2.0*A143.0> 

GO  TO  1520 
1510  21  *  1./2 

XO  «  <  < (-188. 1357*214109. 11 42 1*2 1-23. 79333 >*2142. 050931 >*21 
XO  -  ( ( XO-0. 1 730503 >*2140. 7034845 ) *2 1-0. 0641 09E -3 
XI  -  ( ( (-5.817517*2142. 105874 )*21 -.6896196 )*2l 4 .4952024 X *21 
XI  -  (Xl-O. 187344E-2)*2140. 7979095 

VAL  -  < 1 .O-SQRT (2 ) * ( X1*C0S ( 2-FI/4 . )-X0*SIN<2-FI/4 .))>/( X*X ) 
1520  RETURN 
ENC 
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FORTRAN  IV  V01C-03G  FRI  19-0CT-79  I0JOOSO7 


FACE  OOI 


0001 

0002 

0003 

0005 

0006 

0007 

0008 

0009 

0010  100 
0011 
0012 
0013 

0014 

0015 

0016 

0017 


SUBROUTINE  RFCPS1 
COMMON  /C2/  Z»  R 
IF<Z  .GT.  3)  GO  TO  100 
Z2  «  Z**2 

01  «  (<<2. IE-11  9Z2-5. 38E-9XZ2+6. 757E-7  )*Z2-5. 48443E-5)  *Z2 
01  «  < <01+2.60415E-3>*Z2-6.25E-2>*Z2+.5 
R  -  Z«Q1 
RETURN 
Z1  *  l./Z 

03  *=  <<<-. 14604057*Z1+.27617679)*Z1-. 20210391 )*Z1 +4 . 61835E-3 )*Z1 
03  -  <<03+. 14937 >*Z 1+4. 68E-6  >*Z1  + • 79788456 
04  -  < <<-.21262014*Zl+.19397232)fZl+6.022188E-2)*ZI 
*  -1 . 7222733E-1 >*Z1 

04  -  ( <04+5. 085E-4 ) *Z14 . 37498836 >>Z 1-2. 3561 9449+ Z 
R  *  Q3SC0S ( 04  > /SORT ( Z  > 

RETURN 

ENG 
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FORTRAN  IV 


V01C-O3G  FRl  19-0CT-79  10104:30 


PAGE  001 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0000 

0009 

0010 

0011 

0012 


0013 

0014 

0015 


0016 

0017 


0010 

0019 

0020 

0021 

0022 


0023 

0024 

0025 

0027 

0028 

0029 

0030 

0031 

0032 


0033 


C..  PROGRAM  EDTRFC.FOR 
C.  . 

C..  ALLOWS  INPUT  PARAMETERS  FOR  MULTILAYERED 
C..  REFLECTION  COIL  PROGRAMS  TO  BE  CHANGED 
C. . 

C..  D.  T.  HAYFORD  *  CHIP  WILSON 

C..  BATTELLE  COLUMBUS  LABORATORIES  12-19-78 

C. . 

REAL *4  MATL 
LOGICAL tl  IANS  .  YES 

DIMENSION  COIL ( 13  >  >  AHF'<6).  MATLllO.  3).  ITITLC20) 

DIMENSION  VERB(8).  VERB1 (6) »  VERB2(14>r  VERB3(5). 

•  VERB4  <  7 ) 

DATA  VERB/ 'COIL ' »  'MATL'.  'AMP  '.  'EXIT'.  5*'EXIT'/ 

DATA  VERB1/ 'SAVE'.  ' CHNG ' .  'LIST'.  'HELP'.  'RET  'EXIT'/ 
DATA  VERB2/-R1  '.  'R2  '.  'R3  '»  'R4  '.  'L2  '.  'L3  '. 

*  'L4  ' .  'L5  '.  'L6  '.  'N3  '.  'N4  ',  'RES1'.  'RES2' . 

DATA  VERB3/ ' T  '.  'MU  '»  'RHO'»  'NUM  '.  'NONE'/ 

DATA  VERB4/ 'RES3 ' .  'RES4 ' .  ' CAP3 ' >  'CAP4 ' .  'VOUT',  'GAIN'. 
DATA  YES/'Y'/ 

10  CALL  COMAND ( VERB .  8.  IFILE.  'ENTER  FILE  TYPE?  ') 

GO  TO  (100.  200.  300.  400>  IFILE 

C. . 

C..  'COIL'— INPUT  COIL  PARAMETERS 
C.. 

100  CALL  ATTINC2.  IFILE.  0.  ICUR) 

C..  READ  THE  FILE 

READ < 2 )  ITITL.  COIL 
CALL  CLOSE* 2 > 

C.  . 

C..  GET  NEXT  COMMAND 
110  CALL  COMAND <VERB1.  6.  J.  'COMMAND:  ') 

GO  TO  (120.  130.  140.  150.  160.  400)  J 

C  •  • 

C..  'SAVE'— WRITE  THE  COIL  PARAMETERS  OT  TO  DISK 
120  CALL  ATTOUT < 2.  IFILE.  0.  ICUR> 

CALL  TITLE< ITITL  > 

MRITE<2)  ITITL.  COIL 
CALL  CLOSE  <  2 ) 

GO  TO  110 


C  •  • 

C..  'CHNG'— CHANGE  THE  COIL  PARAMETERS 
130  K  *  1 

CALL  CCMAND ( VERB2 .  14.  Jr  'ITEM!  ') 

IF ( J  . EO .  14)  K  »  2 
GO  TO  (135.  110)  K 

135  TYPE  136.  COIL(J) 

136  FORMAT ( '  OLD  VALUE:  '.  1PE 1 2 . 5/ ' *NEW  VALUE:  ') 
ACCEPT  137.  COIL(J) 

*37  FORMAT (F 10 • 0 ) 

GO  TO  130 


C  <  • 

C..  'LIST'— LIST  PRESENT  COIL  FILE 
140  TYPE  141.  ITITL.  COIL 


'  NONE ' / 
'  NONE '  / 
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PAGC  002 


HI  FORMATUX.  20A2//T10.  'R1  «  '.  F7.4,  T25.  'R2  »  '. 

1  F7.4/T10.  'R3  -  ',  F7 . 4 i  T25.  'R4  «  F7.4/ 

2  T10.  'L2  ■  '.  F7.4»  T25.  'L3  *  '.  F7.4/T10.  'L4  »  ', 

3  F7<  4.  T25.  'L5  «  '•  F7.4/T10.  'L6  «  '.  F7.4/T10.  ' N3 

4  F7  •  1  •  T25.  'HA  *  '»  F7.1/T8.  'RES1  *  F7.1.  T23. 

5  'RES2  -  ' .  F7.1//J 
60  TO  >10 

.  .  'HELP'— GIVE  A  LIST  OF  VARIABLE  NAMES 

150  TYPE  151 

151  FORMAT <  '  Rl  -  DRIVER  COIL  INNER  RAMUS'/ 

•  '  R2  *  DRIVER  COIL  OUTER  RADIUS'/ 

•  '  R3  •=  PICKUP  COIL  INNER  RADIUS'/ 

•  '  R4  «  PICKUP  COIL  OUTER  RADIUS'/ 

•  '  L2  «  LIFTOFF  INCREMENT'/ 

•  '  L3  *  DRIVER  COIL  LENGTH'/ 

•  '  L4  *  PICKUP  COIL  LENGTH'/ 

•  '  L5  -  PICKUP  COIL  RECESS'/ 

•  '  LA  *  MINIMUM  LIFTOFF'/ 

•  '  N3  «  NO.  OF  TURNS-DRIVER' / 

«  '  N4  «  NO.  OF  TURNS-PICKUP'/ 

•  '  RE SI  -  RESISTANCE  OF  DRIVER'/ 

•  '  RES2  -  RESISTANCE  OF  PICKUP'//! 

GO  TO  110 

•  • 

.  .  'RET  '--RETURN  TO  PROGRAM  BEGINNING 
160  GO  TO  10 

.  .  'NATL'  —  INPUT  MATERIAL  PARAMETERS 

I  • 

200  CALL  ATTINC2.  IFILE »  0»  ICUR) 

.  .  READ  THE  FILE 

READ ( 2 )  ITITL »  NUM,  MATL 
CALL  CLOSE ( 2  > 

•  4 

.  .  GET  THE  NEXT  COMMAND 
210  CALL  COMANIK  VERB1 *  6.  J.  'COMMAND!  ') 

GO  TO  <220.  230.  240.  250.  160.  400)  J 

I  • 

,  .  ' SAVE ' --URI TE  MATERIAL  PARAMETERS  OUT  TO  DISK 

220  CALL  ATTOUT  <2.  IFILE.  0.  ICUR> 

CALL  T I TLE (ITITL) 

MRITE<2>  ITITL.  NUM.  MATL 
CALL  CLOSE  <  2 ) 

GO  TO  210 

I  t 

.  .  'CHNG' --CHANGE  MATERIAL  PARAMETERS 

230  CALL  COMANIK  VERB3.  5.  J.  'ITEM!  ') 

GO  TO  <231.  231.  231.  235.  210)  J 

231  TYPE  232 

232  FORMAT! '*L AYER  NO.!  ') 

ACCEPT  233.  K 

233  FORMAT  < 12 ) 

K  -  Ml 
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FORTRAN  IV  V01C-03G  FRI  19-OCT-7?  10:04:30  PAGC  003 

0057  IF(N  .LT .  2  .OR.  N  .GT.  NUH-l)  GO  TO  231 

0059  TYRE  234.  MATL(K»  J> 

0040  234  FORMAT ('  OLD  VALUE:  '*  1  PE  12. 5/ ' *NEU  VALUE:  ') 

0041  ACCEPT  137.  MATL(K*  J) 

0042  GO  TO  230 

0043  235  TYPE  234.  NUM  -  2 

0044  234  FORMAT (  '  OLD  VALUE.’  X2/'*NEW  VALUE:  ') 

0045  ACCEPT  233.  NUMHUM 

0044  NUMHUM  «  NUMHUM  +  2 

0047  IF (NUMHUM  .LT.  1  .OR.  NUMHUM  .GT.  10)  GO  TO  235 

0049  NUM  -  NUMHUM 

0070  IF  (NUM  .ER.  10)  GO  TO  230 

0072  DO  237  I«l»  3 

0073  HO  237  J-NUM+1.  10 

0074  237  MATL< J.  I)  -  0. 

0075  MATL  < 1 .  1)  -  1.E10 

0074  MATL ( 1 •  2)  ■  1. 

0077  MATL < 1 .  3)  «  1.E10 

0076  MATL (NUM.  1 )  *  1 . 

0079  MATL (NUM.  2)  =  1. 

0080  MATL (NUM.  3)  =  1.E10 

0081  GO  TO  230 

C.  . 

C..  'LIST'— LIST  PRESENT  MATERIAL  FILE 
0082  240  TYPE  241,  ITITL.  NUM,  (J,(  MATL(J.  I),  1=1,  3>,  J-l.  10) 

0083  241  FORMAT (IX,  20A2//'  NO.  OF  LAYERS:  '.  12// 

*  T2.  'LAYER'.  T20 .  'THICKNESS'.  T44,  'MU',  T44.  'RHO'// 

*  10(T4.  12.  T20.  1PE10.4,  T40,  E10.4,  T40.  E10.4/>> 

0084  GO  TO  210 

c «  • 

C..  'HELP'— GIVE  LIST  OF  INPUT  PARAMETERS 
0085  250  TYPE  251 

0084  251  FORMAT ( '  LAYERS  ARE  NUMBERED  FROM  DNE  TO  •NUM+2*.  STARTING'/ 

*  '  FROM  THE  BOTTOM.  LAYER  1  IS  ALWAYS  AIR,  WITH  AN  INFINITE'/ 

*  '  THICKNESS.  AND  LAYER'NUM42*  IS  AIR  UI1H  A  THICHNESS  OF  1.'// 

*  'THE  VARIABLES  TO  BE  ENTERED  ARE:'// 

*  '  T  ■  THICKNESS  (IN  INCHES)'/ 

*  'MU  -  RELATIVE  PERMEABILITY'/ 

*  '  RHO  -  RESISTIVITY  (IN  MU-OHM-CM)') 

0087  GO  TO  210 

C.  . 

C..  'AMP  '  —  INPUT  AMPLIFIER  PARAMETERS 
C . . 

0088  300  CALL  ATTIN< 2 »  IFXLE >  0,  ICUR) 

C..  READ  THE  FILE 

0089  READ ( 2 )  ITITL.  AMP 

0090  CALL  CLOSE (2) 

C..  GET  NEXT  COMMAND 

0091  310  CALL  COMAND < VERB 1 ,  4.  J.  'COMMAND:  ') 

0092  GO  TO  (320.  330,  340,  350.  140.  400)  J 

C.  . 

C..  'SAVE'— WRITE  AMPLIFIER  PARAMETERS  OUT  TO  DISK 
0093  320  CALL  ATT0UT(2.  IFXLE,  0.  ICUR) 

0094  CALL  TITLE< ITITL) 
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FORTRAN  IV  V01C-03G  FRI  19-OCT-79  10:04:30  FA6E  004 

0095  WRITE(2)  ITITL.  AMP 

0096  CALL  CLOSE (2) 

0097  GO  TO  310 

C. . 

C..  ' CHNG ' — CHANGE  AMPLIFIER  PARAMETERS 

0098  330  K  *  1 

0099  CALL  C0MAND(VERB4.  7.  Jf  'ITEM:  ') 

0100  IF  <J  .EQ.  7)  K  =  2 

0102  GO  TO  (331.  310)  K 

0103  331  TYPE  234.  AMP(J> 

0104  ACCEPT  137.  AMP(J> 

0105  GO  TO  330 

C. . 

C..  'LIST'— LIST  PRESENT  AMPLIFIER  FILE 
0106  340  TYPE  341.  ITITL.  AMP 

0107  341  FORMAT( IX.  20A2//T10.  'RES3  =  '.  1PE12.5/ 

*  T10.  'RES4  -  '.  E12.5/T10.  'CAP3  *=  '.  E12.5/ 

*  T10.  'CAP4  -  '.  E12.5/T10.  'VOUT  *  '»  E12.5/ 

*  T10.  'GAIN  «  '.  E12.5//J 

0108  GO  TO  310 

C.  . 

C..  'HELP ' --GIVE  LIST  OF  VARIABLE  NAMES 
0109  350  TYPE  351 

0110  351  FORMAT:'  RES3  «  DRIVER  AMPLIFIER  SERIES  RES.  (OHMS)'/ 

*  '  RES4  -  PICKUP  AMPLIFIER  SERIES  RES.  (OHMS)'/ 

*  '  CAP3  -  DRIVER  CIRCUIT  SHUNT  CAP.  (FARADS)'/ 

*  '  CAP4  »  PICKUP  CIRCUIT  SHUNT  CAP.  (FARADS)'/ 

*  '  VOUT  -  DRIVER  OUTPUT  VOLTAGE  (VOLTS)'/ 

*  '  GAIN  «  PICKUP  AMPLIFIER  GAIN'//) 

0111  GO  TO  310 

C.  . 

C..  'EXIT'— EXIT  THE  PROGRAM 

c  • « 

0112  400  TYPE  401 

0113  401  FORMAT < ' *ARE  YOU  SURE  ( Y/N) ?  ') 

0114  ACCEPT  402.  IANS 

0115  402  FORMAT ( A1 ) 

0116  IF( IANS  . EG .  YES)  GO  TO  410 

0118  GO  TO  (110.  210.  310.  10)  IFILE 

0119  410  CALL  EXIT 

0120  END 
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0001 


COO 2 
0003 
0004 
0000 


0006 

0007 

oooe 


ooov 

0010 

0011 

0012 

0013 

0015 

0017 

0018 

0019 

0020 


0021 

0023 

0024 


SUBROUTINE  ATTINO.UN.  NENTRY •  ’FLAG .  I OPEN > 

C.  . 

C..  THIS  SUBROUTINE  ATTACHES  A  f'L £  FOR  INPUT  BY 
C,.  CHECKING  THE  FILE  NAME  DIRECTORY  FOR  THE  FlLENAnL 
C..  AND  DEFAULT  SEQUENCE  NUMBER.  EITHER  THE  DEFAULT 
C. .  MILL  BE  USED  OR  THE  OPERATOR  'ILL  BE  PROMPTED  f OR 

C . ,  ONE . 

C.  . 

C..  INPUT  parameters: 
c. . 

c..  LUN  --  LOGICAL  UNIT  NUMBER  TP  BE  USED  FOR 
C-  THE  FILE  TO  PE  ATTACHED. 

C.  . 

C.  .  NENTRY  --  NUMBER  FOR  THE  PI RECTORY  ENTRY 
C..  TO  BE  OPENED. 

C.  ■ 

C.  .  I  FLAG  —  ~0  PROMPT  OPERATOR  mi*  tf.OUCh'CF.  NO. 

c..  ®-i  USE  DEFAULT  Vol.UE 

C.  •  .GT.  0  USE  IF!. AO  AS  CURRENT  FILE  NUMBER 

C.  . 

C  .  OUTPUT  PARAMETERS 
C  .  . 

r  .  .  ?-Vf.'N  -  GET  EQUAL  W  THE  S’: 0U<: NOE  WO . 

C.  .  OF  1  HE  OPENED . 

C.  . 

LOGICAL*’  STRING! 421 

COMMON/? I LNAh/  LUNUM.  K.  STRING.  I CUR.  I NETT .  1AV 
LUN'JM  *  LUN 
N  =  NENTRY 

C.  . 

C..  ATTACH  THE  FILE  NAME  DIRECTORY . 

C.  . 

10  CALL  DIROPN 
I  *=  I FLAG 

IF(IFLAG)  100.  32.  50 

C.  . 

C.  .  USE  CURRENT  FILE  SEQUENCE  NUMBER*' 

C.  . 

31  IYPL  40;  (STRINGED.  I-'. 2,  14).  ICUK 
40  FORMAT  (' *UrEN  'i  3A:  .  '  FILE  •  13.  '  *»  '  > 

RE AIK  5.  30.  ERR” 35 )  1 
30  FORMAT (110) 

if (i  ,lt,  o  .or.  i  .0".  9??-  nn  tc  35 
IF  (I  .EG.  0)  I  «  I  CUR- 
50  I CUR  *  I 

CALL  DIRCLO 
GO  TO  140 

100  CALL  CLOSE (LUNUM) 

C.  . 

C..  CHECK  FOP  EXISTENCE  OF  T  K 
C  .  . 

140  IF (  IEXIST (STRING)  . GT .  0;  GO  )  ’SO 
TYPE  65.  (STRJNG(I).  1=12.  14).  ICUK 
65  FORMAT ('  3A1.  '  FILE  « ' .  13.  '  NOT  FOUND. ' ) 
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0025 


0026 

0028 

0029 

0030 

0031 

0032 


C. . 

c. . 
c. . 

150 

60 


FRI  19-0CT-79  10J01U2 


V01C-03C 
GO  TO  10 
OPEN  THE  PILE. 

FORMAT?'  *?E4a?>  yVp!,A0'  * STRING < I) »  isl2f 

^PENAfSic^LUN^’  ST^«e.',iJf,'oLrf£1’* '* 

RETURN 

END 


lA)t 
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0001  SUBROUTINE  ATTOUT (LUN.  NENTRY.  IFLAG.  INEW) 

C.  . 

C..  THIS  SUBROUTINE  ATTACHES  A  FILE  FOR  OUTPUT 
C..  BY  CHECKING  THE  FILE  NAME  DIRECTORY  FOR  THE 
C..  FILENAME  AND  THE  NEXT  AVAILABLE  SEQUENCE 
c..  NUMBER;  THEN .  EITHER  the  default  number  will 
C..  BE  USED  OR  THE  OPERATOR  WILL  BE  PROMPTED  FUR 
C..  ONE. 

C.  . 

c..  input: 
c « . 

C..  LUN  —  LOGICAL  UNIT  NUMBER  FOR  THE  FILE 
C..  TO  BE  ATTACHED. 

C. . 

C..  NENTRY  —  THE  DIRECTORY  TO  BE  OPENED. 

C. . 

C..  IFLAG  —  -0  PROMPT  OPERATOR  FOR  SEO.  NO. 

C..  —1  USE  NEXT  AVAILABLE  SEQ.  NO. 

C..  .GT.  0  USE  IFLAG  AS  FILE  SEQ.  NO. 

C..  OUTPUT: 

C.  . 

c..  INEW  --  EQUAL  TO  THE  SEQUENCE  NO.  OF  FILE  OPENED. 

C.. 

0002  L0GICAL*1  STRINC(42> 

0003  COMMON/FILNAM/  LUNUM.  N.  STRING.  I CUR.  I NEXT.  IAV 

0004  LUNUM  -  LUN 

0005  N  -  NENTRY 

C.  . 

C..  ATTACH  THE  FILE  NAME  DIRECTORY. 

C.. 

0006  CALL  DIROPN 

0007  I  -  IFLAG 

0008  IF (IFLAG)  100.  35.  SO 

C.  . 

C..  USE  NEXT  SEQUENCE  NUMBER? 

C.  • 

0009  35  TYPE  40.  (STRING(I).  1-12.  14).  INEXT 

0010  40  FORMAT ( ' tSAVE  AS  '.  3A1 .  '  FILE  ♦  13.  '  ?  ') 

0011  READ (5.  30.  ERR- 35)  I 

0012  30  FORMAT (110) 

0013  IF < I  .LT.  0  .OR.  I  .GT.  999)  GO  TO  35 

0015  50  I CUR  -  I 

0016  IF  <  < I  .NE.  0)  .AND.  (I  .NE.  INEXT))  GO  TO  140 

C .  . 

C..  USE  NEXT  AVAILABLE  SEQUENCE  NUMBER. 

C « » 

0018  100  I CUR  -  INEXT 

0019  INEXT  -  INEXT  ♦  1 

C . . 

C..  CHECK  FOR  EXISTENCE  OF  FILE 
C . . 

0020  140  CALL  DIRCLO 

0021  IFdEXIST(STRING)  .GE.  0)  TYPE  65.  (STRING: I > .  1-12.  14).  ICUR 

0023  65  FORMAT:'  OLD  '.  3A1.  '  FILE  •'»  13.  'DELETED.') 
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f **  M-0C1-79  10101 142 


RETURN 

ENIi 


RRBE  OO; 
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RAGE  OOI 


0001 

C.  . 
C.  . 
C.. 
C.. 
C.. 
C.  . 
C. . 
C.. 

c. . 
c. . 
c.. 
c. . 
c  •  • 
c. . 
c.. 
c. . 
c. . 
c. . 
c  *  • 

0002 

c. . 

0003 

c* « 

0004 

C.. 

c  ♦  • 
c.. 
c.  • 

ooos 

0006 

0007 

Ci  i 

c  • « 
c  •  • 

oooe 

c.. 
c.. 
c. . 

0009 

0010 

001 1  30 

0012 
0013 


SUBROUTINE  PIROPN 

THIS  SUBROUTINE  ATTACHES  THE  FILE  NAME  DIRECTORY, 
LOOKS  UR  THE  REQUESTED  ENTRY,  AND  REPORTS  THE 
CURRENT  AND  NEXT  SEQUENCE  NUMBERS. 

INPUTS 

LUNUM  —  LOGICAL  UNIT  NUMBER. 

NENTRY  —  DIRECTORY  ENTRY  TO  BE  OPENED. 

OUTPUT  S 

STRING  —  ARRAY  CONTAINING  THE  FILE  NAME. 

I CUR  —  CURRENT  SEQUENCE  NUMBER. 

INEXT  —  NEXT  SEQUENCE  NUMBER. 

L0GICALB1  STRING (42) ,  DUM1(3),  DUM2(3> 

COMMON/FILNAM/  LUNUM,  NENTRY,  STRING,  ICUR,  INEXT,  IAV 
EQUIVALENCE  (DUMl(l),  STRING(S)),  (DUM2<1),  STRINGS  16) > 


ATTACH  THE  FILE  NAME  DIRECTORY. 

CALL  ASSIGNiLUNUM,  'DIRS RFCRFC . DIR ' ,  14) 
DEFINE  FILE  LUNUM<4,  21,  U,  IAV) 

IAV  -  NENTRY  +  1 

READ  IN  THE  REQUESTED  ENTRY. 

READ(LUNUM'IAV)  STRING 

DETERMINE  CURRENT  AND  NEXT  AVAILABLE  SEQ.  NO. 

DECODE  <3,  30,  DUM1 )  ICUR 
DECODE  (3,  30,  DUM2)  INEXT 
FORMAT  < 1 3 ) 

RETURN 

END 
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FORTRAN  IV  V01C-03G  FRI  19-0CT-7?  10102130  FACE  001 

0001  SUBROUTINE  DIRCLO 

Cat 

Ca.  THIS  SUBROUTINE  UPDATES  AND  CLOSES  THE  FILE 
C..  HAHE  DIRECTORY . 

Ca  a 

Ca.  INPUT: 

C.  a 

Ca.  LUNUM  —  LOGICAL  UNIT  NUMBER. 

C.  a 

Ca.  NENTRY  —  DIRECTORY  ENTRY  TO  BE  UPDATED. 

Ca. 

Ca.  STRING  —  ARRAY  CONTAINING  THE  FILENAME. 

C.« 

Ca.  I CUR  —  CURRENT  SEQUENCE  NUMBER. 

C  a  a 

0002  L0GICALB1  STRING<42>,  DUM1<3)»  DUM2<3>»  DUM3 >  DUMA 

Ca. 

0003  COMMON/FILNAM/  LUNUM.  NENTRY ,  STRING.  ICUR.  I NEXT .  IAV 

Ca. 

0004  EQUIVALENCE  <DUMlU)a  STRING<B>>»  (DUM2<1>.  STRING<16>> 

Ca. 

C.  . 

Ca.  ENCODE  ICUR  AND  I NEXT  INTO  STRING 

C  a  a 
C  a  a 

0005  DUM3  *  STRING ( 11) 

0006  ENCODE  (3.  30.  DUM1 )  ICUR 

0007  STRING(ll)  -  DUM3 

C  a  a 

0008  DUM3  -  STRING< 19) 

000?  ENCODE  <3.  30,  DUM2)  I NEXT 

0010  STRINGU?)  -  DUM3 

C.  . 

0011  30  FORMAT (I 3) 

C  a  a 

C.a  CONVERT  BLANKS  TO  ZEROS. 

C  a  a 


0012 

IF (STRING(6 ) 

.EQ. 

*040) 

STRING(B) 

K 

*060 

0014 

IF  <STRING<9) 

.EQ. 

*040) 

STRING(?) 

m 

*060 

0016 

IF (STRING( 16) 

.EQ. 

*040) 

STRING! 16) 

m 

*060 

0018 

IF ( STRING< 17  > 

.EQ. 

*040) 

STRING < 17) 

m 

•060 

0020  IAV  «  NENTRY  ♦  1 

C  a  a 

Ca.  UPDATE  DIRECTORY  ENTRY  AND  CLOSE  IT. 

C  a  a 

0021  WRITE <LUNUM' IAV)  STRING 

0022  CALL  CLOSE (LUNUM) 

0023  RETURN 

0024  END 
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0001 


0002 

0003 

0004 

0005 

0006 

0007 

0009 

0010 

0011 

0012 

0013 

0014 
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RACE  001 


C  4  • 

c. . 
c. . 
c. . 
c. . 
c. . 
c. . 
c  • « 


SUBROUTINE  TITLE!  ITITLE) 
PROGRAMMED  BY  CHIP  WILSON. 


THIS  SUBROUTINE  LISTS  THE  40  CHARACTER  TITLE  CONTAINED  IN 
ARRAY  ITITLE f  ASKS  THE  OPERATOR  IF  HE  HANTS  TO  CHANGE  ITf 
AND  RETURNS  WITH  EITHER  THE  OLD  TITLE  OR  A  NEW  TITLE  AS 
ENTERED  BY  THE  OPERATOR. 


LOGICALtl  ITITLE (40) f  ANS 

c«  • 


TYPE  200 f  ITITLE 

200  FORMAT! '  PRESENT  TITLE!  >'f  40A1f  'C'f  /f  '•NEW  TITLE  (Y/N)? 

ACCEPT  21 Of  ANS 
210  FORMAT!  40A1 ) 

IF!  ANS  .NE.  *131)  RETURN 
DO  300  J*1 f  40 
300  ITITLE ! J)  •  *040 

CALL  MSGOUT! 'ENTER  NEW  TITLES  ') 

ACCEPT  210f  ITITLE 

RETURN 

END 


'  > 
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FORTRAN  IV  V01C-03G  FRI  19-0CT-79  10:03:08  PAGE  001 

0001  SUBROUTINE  COMAND  <  ARRAY  *  ISIZE.  HATCH.  PROMPT) 

C  •  • 

C..  PROGRAMMED  BY  CHIP  WILSON. 

C  •  .  * 

C..  A  FOUR  CHARACTER  COMMAND  IS  INPUT  FROM  THE 
C..  KEYBOARD  AND  A  8EARCH  IS  MADE  FOR  A  MATCH. 

C..  A  NONAMBIGUOUS  ABREVIAT10N  IS  CONSIDERED  A 
C..  VALID  MATCH. 

C..  IF  NO  MATCH  IS  FOUND.  ALL  COMMANDS  ARC  LISTED  OUT 
C..  AND  THE  ROUTINE  WAITS  FOR  ANOTHER  INPUT. 

C. . 

C..  INPUT  - 

C..  ARRAY  —  ARRAY  OF  COMMANDS  TO  BE  CHECKED. 

C..  ISIZE  —  NUMBER  OF  COMMANDS  TO  BE  CHECKED. 

C..  PROMPT  —  CHARACTER  STRING  CONTAINING  A  PROMPTING 
C..  MESSAGE  TERMINATED  BY  A  NULL. 

C . . 

C..  OUTPUT  - 

C..  HATCH  —  ARRAY  SUBSCRIPT  WHERE  HATCH  OCCURED. 

C.  . 

0002  L0GICALF1  PROMPT  Cl).  LTRYC4) .  REPLYC4) 

0003  EQUIVALENCE  (TRY.  LTRYC 1 ) > .  (ANSWER.  REPLY(l)) 

0004  DIMENSION  ARRAY ( 1 ) 

0005  3  DO  50  J-1.40 

0006  IF (PROMPT ( J)  .EG.  0)  GO  TO  5 

0008  50  CONTINUE 

0009  5  TYPE  4.  (PROMPT ( I >. I»1 . J> 

0010  4  FORMAT ( /.  40A1) 

0011  READ ( 5 >  110.  ERR-3)  ANSWER 

0012  110  FORMAT (A4> 

0013  IFIND  ■  0 

0014  LETRS  -  0 

0015  DO  10  K-l.  ISIZE 

0016  TRY  ■  ARRAY (K ) 

0017  DO  9  J«l.  4 

0018  IF ( REPLY ( J)  .NE.  LTRY(J))  GO  TO  20 

0020  9  CONTINUE 

0021  20  IF( J-l  -  LETRS)  10.22.21 

0022  21  IFIND  -  0 

0023  LETRS  -  J-l 

0024  22  IFIND  »  IFIND  +  1 

0025  MATCH  ■  K 

0026  10  CONTINUE 

0027  IF( IFIND. EG. 1. AND. (LETRS. EG. 4  .OR.  REPLY (LETRS41 ) .EG. *040) >  RETURN 

0029  TYPE  210.  ( ARRAY ( I ) .  I  -  1.  ISIZE) 

0030  210  FORMAT  ('♦'»  4(  10(A4J .  ')/.  '  ')> 

0031  GO  TO  3 

0032  END 
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FACE  001 


0001 


0002 

0003 


0004 

0005 

0006 

0007 

0008 

0009 

0011 

0012 

0013 

0014 

0015 


C.  . 

c«  • 
c. . 
c. . 
c. . 
c. . 
c. . 
c. . 
c. . 
c. . 
c. . 


c. . 
c. . 


FUNCTION  IEXIST (FNAME) 

THIS  FUNCTION  PERFORMS  A  DIRECTORY  LOOK  UP 
TO  CHECK  FOR  THE  EXISTENCE  OF  A  FILE. 

IEXIST  -  -1  ERROR. 

■  -2  FILE  NOT  FOUND. 

.GE.  0  FILE  FOUND.  IEXIST  •  NO.  OF  BLOCKS  IN  FILE. 

FNAHE  ■  A  14-CHARACTER  LOGICAL  ARRAY  CONTAINING  THE  FILE  NAME 
IN  ASCII.  E.G.  DEV1FLNAME.EXT 

INTEGER  RAD50 ( 4 ) 

L0GICALS1  FNAME(l) 

CONVERT  FILENAME  TO  RADIX  50  FORMAT 
CALL  IRAD50< 3*  FNAME<1).  RAD50(1>> 

CALL  IRAD50<6>  FNAME < 5 > »  RAD50(2>> 

CALL  IRAD50(3f  FNAME < 12) *  RADS0(4>> 

IEXIST  -  -1 

1CHAN  -  IGETCO 

IF< I CHAN  ,LT .  0)  RETURN 

IEXIST  -  LOOKUP: I CHAN.  RAD50 > 

CALL  CLOSEC(ICHAN) 

CALL  IFREEC(ICHAN) 

RETURN 

END 
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0001  SUBROUTINE  HSGOUT(NSG) 

C.. 

C..  PROGRAMMED  BT  CHIR  WILSON 
C.. 

C..  THIS  ROUTINE  SENDS  AN  ASCII  ME6SAGE  OUT  TO  THE 

C..  TERMINAL.  THE  L0GICAL91  STRING*  MSG*  MUST  BE 
C..  TERMINATED  BT  A  NULL  BTTE. 

C.. 

0002  LOGICALtl  MSG(l) 

C. . 

0003  BO  10  J-l.  SO 

0004  IF  <  MSG(J)  .EO.  0  >  RETURN 

0006  S  IF  (  ITT  OUR  <  MSG ( J ) )  .EO.  1  )  GO  TO  5 

0008  10  CONTINUE 

0009  RETURN 

0010  END 
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